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Engineering  Division,  EL.  and  Dr.  John  Harrison,  Director.  EL. 

Dr.  John  O.  Curtis  collected  the  information  and  prepared  the  report  under 
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nical  Report  EL-93-25,  U.S.  Army  Engineer  Waterways  Experi¬ 
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Conversion  Factors,  Non-SI  to 
SI  Units  of  Measurement 


Non-Sl  units  of  measurement  used  in  this  report  can  be  converted  to  SI  units 
as  follows: 


Multiply 

By 

To  Obtain 

degrees  (angle) 

0  01745329 

radians 

feet 

0  3046 

meters 

inches 

254 

cenDmeters 

1  introduction 


New  and  improved  methods  of  remote  sensing  have  increased  the  under¬ 
standing  of  Earth's  origins,  its  resources,  and  those  processes  that  contribute  to 
its  dynamic  (on  a  large  time  scale)  nature.  Scientists  and  engineers  from  many 
discipLnes  are  con.stanlly  exploring  new  ways  to  quantity  the  Earth’s  properties 
for  ihe’r  particular  applications.  Among  these  methods  arc  measurements  of 
electromagnetic  energy  in  many  different  wavelength  regimes,  both  passive 
and  active  The  development  of  small  powerful  sources  and  ultrasensitive 
receivers  along  with  improved  data  processing  capabilities  has  fostered 
renewed  interest  in  measurements  of  the  Earth  within  the  microwave  region  of 
the  spectmm,  where  wavelengths  in  air  range  from  a  few  millimeters  to  several 
meters. 

Natural  terrain  surfaces  consist  of  bare  soils,  rocks,  vegetation,  atid  water 
This  study  was  initiated  with  the  hope  of  making  a  meaningful  contribution  to 
the  understanding  of  microwave  interactions  with  natural  terrain.  In  particular, 
the  focus  here  is  to  summarize  what  is  currently  known  about  the  measurement 
and  modeling  of  the  elccuncal  properties  of  well-characterized  soils.  Follow-on 
reports  will  deal  with  a  new  measurement  capability  recently  developed  at  the 
L'.S.  Army  Engineer  Waterways  Experiment  Station  (N^'ES)  in  Vicksburg,  MS, 
and  the  interpretation  of  data  collected  at  WES  using  those  new  measurement 
resources.  Most  of  the  information  contained  in  this  series  of  reports  is  drawn 
from  a  doctoral  research  program  recently  completed  by  the  author  (Curtis 
1992). 

Reasons  for  Studying  Electrical  Properties  of 
Soils 

Soil  moisture  measurements 

Tlie  dominant  factor  that  controls  the  electrical  behavior  of  soils  is  the  pres¬ 
ence  of  water  (Topp,  Davis,  and  Annan  1980).  Obviously  one  would  hope  to 
take  advantage  of  this  experimental  fact  to  develop  a  means  of  accurately  and 
quickly  measuring  the  moisture  in  soils  without  having  to  collect  numerous 
samples  in  the  field,  weigh  them,  dry  them  for  extended  jseriods  of  time,  and 
weigh  them  again  to  obtain  either  a  gravimetric  (weight  of  water/weighi  of  dry 
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soil)  moisture  content  or  a  more  uscl'ul  volumetric  (volume  of  waicrAolume  of 
soil  sample)  moistua’  content. 

Numerous  attempts  have  been  made  to  develop  a  uselul  method  lor  measur¬ 
ing  soil  moisture  content.  All  have  met  with  varj'ing  degrees  of  success  and 
none  ha.s  proven  accurate  under  all  conditions.  For  example,  a  technique  tliai 
involves  burxing  radar  transmit  and  receive  antennas  in  the  soil  and  relating 
tlic  measurements  of  attenuated  received  signals  to  moistum  content  (Birchak 
ct  al.  1974)  is  destructive  to  the  soil  fabric  (the  way  in  which  soil  panicles  arc 
arranged),  is  very  much  controlled  by  the  fabric  and  the  si/.c  distribution  of 
panicles,  and  precludes  the  u.sc  of  the  same  instrument  in  multiple  locations 
'nature  is  not  homogeneous)  as  well  as  the  ability  to  easily  repair  defective 
equipment. 

Another  approach  for  making  field  measurements  of  soil  moi.sture  taken  by 
some  researchers  (that  also  has  numerous  application  in  the  biomedical  field)  is 
that  of  mcasunng  the  change  in  fnngc  capacitance  of  an  open-ended  coaxial 
probe  (Thomas  1966,  Brunfcldt  1987;  Gabncl,  Grant,  and  Young  1986). 

When  pressed  against  a  soil  whose  propcnics  arc  unknowm,  the  resulting 
change  in  capacitance  pioduccd  by  the  impedance  mismatch  is  related  to  elec¬ 
trical  properties  through  calibration  relationships.  Problems  arising  from  the.sc 
mcasurcmcnis  include  the  need  to  have  proper  contact  between  the  probe  lip 
and  the  soil  surface,  the  fact  that  the  volume  of  material  associated  with  the 
fringe  capacitance  is  quite  small  ( 1  cm’  or  less),  and  that  calibration  conditions 
simply  cannot  account  for  all  of  the  diektiric  loss  mechanisms  that  exist  in 
natural  soils.  The  losses  in  moist  soils  can  be  highly  frequency  dependent 
over  a  range  of  several  frequency  decades  on  the  electromagnetic  spectrum.  A 
recent  variation  on  the  open-ended  probe  measurement  scheme  involves  the 
use  of  a  waveguide  section  instead  of  a  coaxial  device  (Parchomchuk,  WaJlcn- 
der,  and  King  1990). 

The  concern  over  small  sample  volumes  can  be  overcome  with  a  redesign 
of  the  open-ended  coaxial  probe  that  replaces  the  sobd  outer  conductor  with 
several  pointed  lines  (Campbell  1988)  which  allows  for  the  probe  to  be  pushed 
into  the  surface  of  soft  soils.  The  volume  of  soil  enclosed  by  such  probes  can 
easily  be  lens  of  cubic  centimeters.  The  lined  coaxial  probes  that  have  been 
built  ID  date  operate  in  a  frequency  range  that  is  very'  much  subject  to  the 
material-dependent  loss  mechanisms  alluded  to  above.  Rc.scarch  continues  that 
is  directed  to  better  undersunding  these  various  mechanisms  and  to  fabricate  a 
probe  that  operates  in  a  frequency  range  that  is  not  subject  to  such  material - 
dependent  anomalies. 


Subterranean  Investigations 

Many  types  of  electrical  measuromenus  in  soils  (and  rocks)  arc  used  to 
undcr.vtand  what  lies  beneath  the  Earth's  .surface  (Telford  cl  al.  1984).  Among 
these  arc  vanous  schemes  for  measuring  soil  resistivity  a.s  well  as  the 
attenuation  of  propagating  clcciromagnciic  waves. 
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Resistivity  measurements.  Resistivity  data  in  soils  sue  collected  by  inject 
ing  known  currents  (usually  of  a  frequency  less  than  60  Hz;  altcmating  current 
(AC)  required  to  minimize  effects  of  charge  buildup  on  the  probe)  into  the 
ground  and  measuring  potential  differences  across  pairs  of  nearby  electrodes. 
Assuming  homogeneous  media  and  uniform  resistivity,  it  is  possible  to  calcu¬ 
late  from  the  potential  differences  an  apparent  resistivity  of  the  earth  material. 
These  measurements  will  be  affected  by  the  presence  of  water  in  the  soil  or 
rock,  by  the  presence  of  mineral  compounds  that  could  go  into  ionic  solution 
with  available  water,  and  by  the  physical  structure  of  the  subsurface  terrain 
itself.  This  being  the  case,  resistivity  measurements  are  useful  as  a  measure  of 
subsurface  water  volume,  the  locations  of  mineral  deposiLs,  and  subsurface 
structure. 

Electromagnetic  wave  propagation.  Another  method  of  making  electrical 
subsurface  measurements  involves  the  transmission  '^f  electromagnetic  waves 
into  the  soil  and  detection  of  energy  that  results  from  waves  rcllccted  from 
subsurface  anomalies  Having  the  ability  to  track  wave  propagation  in  time, 
either  by  pulsing  the  source  or  sweeping  over  a  known  frequency  band  in 
some  controlled  manner,  means  that  electromagnetic  wave  propagation  meth¬ 
ods  of  subterranean  investigations  are  particularly  useful  for  locating  the  depth 
of  electrical  anomalies  such  as  die  water  table  in  sandy  soil  (Olhoeft  1983; 
Stewart  1982;  Wright.  Olhoeft,  and  Watts  1984),  buried  pipes  or  wires,  or 
cavities  such  as  tunnels  or  caves  (Ballard  1983).  Other  applications  include 
the  delineation  of  stratified  media  (Lundien  1972)  and  determination  of  the 
thickness  of  ice  and  frost  layers  (Jakkula,  ^linen,  and  Tiuri  1980;  O’Neill  and 
Arcone  1991). 

There  are  some  practical  bounds  on  the  utility  of  radio  frequency  systems 
to  conduci  subterranean  investigations  due  to  the  pheromenon  of  "skin  depth" 
(Appendix  .A),  a  measure  of  the  attenuation  of  the  electromagnetic  energy  as  it 
travels  through  the  medium.  For  relatively  low-frequency  sources  (about 
200  MHz),  it  can  be  shown  that  low-loss  soils  such  as  dry  sands  can  possess  a 
skin  depth  of  about  10  to  15  m,  while  high-loss  soils  such  as  wet  silts  and 
clays  may  have  skin  depths  of  only  a  few  cenbmeters.  Modem  radio 
frequency  receivers  are  extremely  sensitive  devices,  often  having  a  dynamic 
range  of  50  to  100  db.  The  signal  at  skin  depth  represents  about  an  8.7-db 
loss  in  power  or  a  two-way  loss  at  the  receiver  for  reflected  signals  of  about 
17.4-db.  It  is  certainly  not  inconceivable  that  radio  frequency  receivers  should 
be  capable  of  successfuUy  delecting  reflected  signals  from  subterranean 
anomalies  at  depths  of  two  to  three  skin  depths  or  more. 


Remote  sensing  of  environment 

Virtually  al!  remote  sensing  of  the  environment  from  airborne  or 
spacebome  platforms  involves  the  measurement  of  electromagnetic  radiation 
from  the  Earth’s  surface  and/or  atmosphere.  Passive  surveillance  involves 
measurements  of  emitted  radiation  and  that  reflected  from  natural  sources  such 
as  the  sun,  the  atmosphere,  and  surrounding  terrain.  Active  remote  .sensing 
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measurement  systems  include  a  source  for  illuminating  the  target  of  interest. 
Whether  passive  or  active,  whether  visual,  thermal  infrared,  microwave,  or 
millimeter  wave  sensors  aic  utilized,  remote  sensing  is  the  collection  and  inter¬ 
pretation  of  electromagnetic  radiation  and,  as  such,  demands  an  understanding 
of  the  dielectric  properties  of  those  materials  being  observed. 

Environmental  remote  sensing  applications  form  a  list  that  grows  yearly  as 
electronic  components  are  improved  and  data  collection  and  processing  hard¬ 
ware  and  software  become  faster,  more  reliable,  and  less  costly.  For  example, 
satellites  can  provide  worldwide  surveys  of  land-use  patterns  to  monitor  the 
threat  of  urtaniza'.ion.  waste  disposal,  and  erosion  of  the  land  (Colwell  1983). 
Similar  systems  (including  those  mounted  in  aircraft)  can  monitor  the  health  of 
vegetation  to  keep  abreast  of  such  things  as  loss  of  forest  and  the  potential  for 
food  shonages.  Sea  traffic  in  the  far  northern  and  southern  .shipping  lanes  can 
be  made  safer  through  the  use  of  airborne  and  spacebome  sensors  to  detect  ice 
hazards. 

One  of  the  more  obvious  applications  for  microwave  remote  .sensing 
devices  is  that  of  conducting  surface  moisture  surveys  to  help  predict  ground- 
water  availability  and  the  potential  for  flooding.  Attempts  have  been  made  to 
relate  soil  moisture  to  both  laboratory  reflectance  data  (Lundien  1966)  and  to 
airborne  sensor  backscatier  measurements  (Ulaby,  Cihlar,  and  Moore  1974; 

John  1992).  Careful  airborne  sensor  measurements  might  provide  a  first 
approximation  to  the  complex  dielectric  constant  of  the  soil  near  the  surface. 

Because  of  the  difference  in  dielectric  behavior  of  liquid  water  and  various 
forms  of  ice,  it  may  be  possible  to  use  airborne  sensors  to  detect  freezing  and 
thawing  (WcgmuUer  1990)  in  remote  locations  that  could  be  used  to  predict 
spring  runoff  conditions  and  all  that  encompasses  for  agriculniraJ  applications, 
the  effects  on  the  fishing  industry,  and  the  anticipation  of  flooding  in  built-up 
areas.  Military  analysts  are  concerned  about  soil  moisture  conditions  because 
of  its  impact  on  irafficability,  the  ability  of  vehicles  to  move  effectively  over 
natural  terrain. 

Another  remote  sensing  application  is  the  mapping  of  exposed  soils  and 
rocks  in  remote  areas  of  the  world  from  high-flying  aircraft  or  satellites,  which 
might  prove  useful  for  geomorphological  studies  (Swanson  et  al.  1988)  or  even 
mineral  exploration.  A  less  obvious  but  recent  application  of  microwave 
remote  sensing  in  soils  dealt  with  archeological  surveys  in  desert  areas  (Berlin 
et  al.  1986;  McCauley  et  al.  1986). 


Others 

While  the  above  paragraphs  emphasize  some  of  the  most  obvious  and  use¬ 
ful  applications  of  a  better  understanding  of  soil  electrical  properiics.  others 
have  been  noted  in  the  literature.  For  example,  some  researchers  have 
attempted  to  relate  electrical  property  measurements  to  the  physical  properties 
of  soils  (Campbell  and  Ulrichs  1969,  Hayrc  1970;  Amlanandan  and  Smith 
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1973;  Madden  1974).  Of  course,  nothing  has  been  said  about  the  military’s 
need  to  better  understand  the  microwave  response  of  soils  that  form  the  back¬ 
grounds  to  military  targets;  i.e.,  when  and  why  does  clutter  become  a  source  of 
target-like  signatures.  Outside  the  topic  of  soils,  studies  of  the  microwave 
response  of  foodstuffs  have  direct  application  to  quality  control  concerns  in  the 
food  industry  (Nelson  1973,  1983). 

Another  new  application  of  this  technology  that  is  closely  related  to  the 
discussion  on  soil  moisture  is  that  of  detecting  liquid  ground  contaminants, 
either  near  the  surface  or  at  arbitrary  depths  using  a  specially  fabricated  probe. 
If,  as  will  be  argued  later,  polarizable  liquids  can  be  characterized  by  a  unique 
frequency  of  peak  losses  due  to  the  dielectric  relaxation  phenomenon,  then  a 
probe  could  be  designed  to  measure  losses  over  a  frequency  span  tiiat  is  broad 
enough  to  detect  a  peak  loss  frequency  and,  coupled  with  the  results  of  a 
thorough  experimental  program,  to  identify  the  particular  contaminant. 


Complex  Dielectric  Constant 

Background  information  for  this  study  would  not  be  complete  without  a 
definition  of  terms  that  will  be  used.  Consider,  for  the  moment.  Ampere's  law 
written  for  linear,  isotropic  materials  and  for  current  density  divided  into  a 
component  due  to  the  motion  of  free  charges  such  as  electrons  and  ions  and 
one  due  to  other  factors  such  as  time  variation  of  polarization  and  magnetiza¬ 
tion  (Gaussian  units  arc  used  throughout  the  text.  For  other  units  see  the 
excellent  Appendix  on  Units  and  Dimensions  in  Jackson  (1975)). 


Vx/7  =  Ill 


_  (^4rc'' 

J  + 

f'1 

dD 

IT 


(1) 


where 

//  =  the  magnetic  field 
c  =  the  speed  of  hghi  in  a  vacuum 
1 1  =  the  current  density  due  to  free  charges 
D  =  the  electric  displacement. 

If  one  also  assumes  that  the  magnetic  con.nbution  to  free  charge  forces  is 
small  relative  to  that  of  the  electric  field,  il.  then  one  can  rewrite  Equation  1 
using  Ohm's  law. 

=  o£  (2) 
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where  o  is  the  conductivity  of  the  medium.  Furthermore,  having  assumed  a 
linear  isotropic  material,  one  can  write  the  electric  displacement  in  terms  of  the 
electric  field  as 

D  ^  e  E  (3) 


where  e  is  the  electric  permittivity  of  the  material.  (It  is  understood  that  if  £  is 
frequency  dependent,  ^uation  3  is  not  rigorously  correct  in  the  time  domain, 
but  rather  that  the  electric  displacement  and  electric  field  are  related  through 
Fourier  transforms  of  the  frequency-dependent  perminivity  (Jackson  1975, 
p.  307).  Adding  this  noutional  complexity  would  only  serve  to  cloud  the 
qualitative  development  intended  within  this  section.)  Then  Ampere’s  law 
can  be  written  as 


V  X  H 


Ana  p  e  3£ 
- £  ♦  — _ 

c  c  dt 


(4) 


Thus,  the  flow  of  current  in  a  medium  is  both  proportional  to  the  applied  elec¬ 
tric  field  and  to  the  time  rate  of  change  of  that  field.  By  analogy  with  simple 
electrical  circuits,  permittivity  is  a  measure  of  the  electrical  capacity  or  capaci¬ 
tance  of  the  media,  that  is,  its  ability  to  store  charge.  Capacitor  plates  in  a 
vacuum  collect  and  give  up  charge;  thus,  there  is  a  permittivity  of  free  space, 

£„  (unity  in  Gaussian  units),  Capacitor  plates  filled  witli  a  polarizable  material 
exhibit  a  higher  capacitance  due  to  material  polarization  such  that  the  material 
permittivity  can  be  written  as 

£  =  e^^l  +  AnX,)  (5) 


where  X,  is  the  electric  susceptibility  for  linear  materials  and  relates  Lhe 
materia]  polarization  to  electric  field.  The  value  of  material  permittivity  nor¬ 
malized  to  the  pennitiivity  of  free  space  is  what  is  called  the  dielectric  con¬ 
stant  (or  relative  permittivity). 

e  s  -i  =  1  ♦  47tX  (6) 

Ed 


From  this  point  on  in  the  text,  e  will  refer  to  the  relative  permittivity. 

As  will  be  shown  by  the  data,  the  diclccuic  constant  of  real  matcriais  is 
frequency  dependent  (a  characteristic  referred  to  as  "dispersion”).  At  relatively 
low  frequencies,  polarizable  panicles,  which  are  effectively  electric  dipoles,  are 
able  to  rotate  or  align  themselves  with  the  changing  electric  field  with  linJe 
loss  in  energy.  However,  as  will  be  discussed  later,  several  other  mechanisms 
arc  at  work  that  do  result  in  electrical  energy  losses,  particularly  at  frequencies 
less  than  1(X)  MHz.  As  the  frequency  goes  up,  the  dipoles  cannot  keep  up 
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completely  with  the  electric  field,  resulting  in  a  phase  shift  in  current  that 
peaks  in  the  10  to  20  GHz  frequency  range  and  is  quite  temperature  depen¬ 
dent.  Finally,  at  frequencies  well  above  the  peak  loss  range,  the  dipoles  do  not 
even  respond,  resulting  in  a  reduction  of  real  capacitance  to  its  optical  value 
that  is  due  to  electronic  polarizability  (FeyTiman,  Leighton,  and  Sands  1964). 
This  frequency-dependent  loss  mechanism  can  be  modeled  by  a  complex 
dielectric  constant 

E  =  e'  iz'' 


where  e"  reflects  the  phase  lag  in  the  particle  motion, 
law  could  be  written  as 

V  ^  u  dE 

y  X  H  ^ -  E  +  _  +  /  - ; 

c  c  c  j  at 


Finally,  then.  Ampere’s 

(8) 


If  one  further  assumes  that  the  material  is  subjected  to  a  sinusoidal  electric 
field 


E  =  E^e 


(9) 


then 


I 


V 


X  //  = 


47tO 

C 


(10) 


The  first  two  terms  on  the  right  hand  side  of  Equation  10  represent  the  con¬ 
ductive  nature  of  the  material,  the  flow  of  energy  that  is  associated  with  losses. 
The  third  term  is  a  measure  of  its  energy  flow  due  to  its  polarizability.  All  of 
the  terms  within  the  square  brackets  are  collectively  referred  to  as  the  admit¬ 
tance  of  the  material,  being  the  ratio  of  current  to  voltage.  The  first  two  terms 
represent  the  conductance  of  the  material,  and  the  last  term  represents  its 
susceptance.  It  is  clear,  then,  that  u£"/4n  can  be  thought  of  as  a  dielectric 
conductivity. 


It  appears  that  electrical  property  measurements  reported  in  the  literature 
seldom,  if  at  all.  distinguish  between  the  two  physically  different  loss  mecha¬ 
nisms.  In  fact,  the  terms  e"  and  4Tto,/co  are  often  used  interchangeably.  This 
observation,  coupled  with  the  fact  that  the  engineering  community  uses  a 
noiational  and  sign  convention  for  dealing  with  dielectric  properties  that  is 
different  from  that  used  by  physicists  and  chemists,  makes  reading  of  the  liter¬ 
ature  often  quite  burdensome. 
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2  Electrical  Properties  of 
Water 


Certainly,  over  the  range  of  the  cleciromagnetic  spectrum  from  direct  cur¬ 
rent  (DC)  to  millimeter  waves,  the  single  most  important  factor  in  determining 
the  electrical  behavior  of  soils  is  the  presence  of  water.  Whether  as  a  solvent 
that  provides  a  medium  and  path  for  ionic  conduction,  or  as  panicles  that  bond 
to  clay  changing  the  electrochemical  characteristics  of  the  soil,  or  as  a  source 
of  dipoles  that  lead  to  high-frequency  residual  dielectric  conduction,  water  has 
a  major  impact  on  how  the  moist  soil  interacts  with  electromagnetic  fields. 

For  these  reasons,  a  study  of  the  electrical  properties  of  moist  soils  should 
include  a  study  of  the  electrical  properties  of  water. 


Summary  of  Data  on  Pure  Water 

The  available  literature  on  electrical  measurements  of  water  show  a 
consensus  on  several  points.  First,  the  so-called  static  (frequency  <1(X)  MHz) 
dielectric  behavior  of  pure,  deionized  water  is  very  well  defined  over  a  broad 
range  of  temperatures  (Hasted  1973).  Second,  there  is  a  loss  mechanism 
(dipole  relaxation)  at  higher  frequencies  (greater  than  1  GHz)  that  can  be 
modeled  in  a  rather  simple  fashion  and  for  which  data  is  highly  repeatable 
(Kjatze  1986).  Third,  fresh  water  exhibits  low-frequency  conductivity  losses. 
Finally,  seawater  and  saline  solutions  have  low- frequency  losses  that  are  a 
couple  of  orders  of  magiutude  higher  than  fresh  water.  TTierc  seems  to  be 
little,  if  any.  information  on  the  electrical  behavior  of  pure,  deionized  water  in 
the  frequency  range  between  100  MHz  and  1  GHz. 

Figure  1  is  a  graphical  representation  of  the  complex  dielectric  constant  for 
liquid  fresh  water  as  a  function  of  temperature.  These  curves  were  plotted 
from  empirical  relationships  developed  by  Ray  (1972).  for  which  excellent 
correlation  between  tlic  empirical  fits  and  actual  data  was  clearly  demoastrated. 
A  description  of  the  models  for  losses  at  higher  frequencies  follows. 
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Figure  1 .  Complex  dielectric  constant  tor  liquid  water  as  a  function  ot  tem¬ 
perature  (analytical  models  after  Ray  (1972)) 


Microwave  Frequency  Loss  Mechanisms 

Breaking  of  hydrogen  bonds 

The  precise  physical  structure  of  liquid  water  is  unknown.  A  physical 
model  of  free  water  molecules  (H2O)  is  accepted  as  consisting  of  two  hydrogen 
atoms  bonded  to  an  oxygen  atom  in  a  V-shape  with  an  angle  of  104.5  deg‘ 
formed  by  lines  through  the  hydrogen  protons  and  meeting  at  the  center  of  the 
oxygen  nucleus.  Viewing  such  a  molecule  along  the  normal  to  the  plane 
through  the  bond  angle  using  the  point  charge  model  of  the  water  molecule 


'  A  table  of  factors  for  converting  non-Sl  units  of  measurement  to  SI  units  is  presented  on 
page  viu. 
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shown  in  Figure  2,  one  can  compute  that  the  free  water  molecule  has  a  dipole 
moment  with  a  value  of  1.83  x  10  '*  electrostatic  units.  The  molecule  has  two 
positively  charged  "comers’’  at  the  hydrogen  protons  and  two  negatively 
charged  comers  at  the  electron  pairs  that  arc  not  shared. 


Figure  2.  Point  charge  model  ol  free  water  molecule  (from  Mitchell  (1976)) 


Hydrogen  bonds  are  formed  when  one  of  these  positive  comers  bonds  to 
the  negative  comer  of  a  neighboring  molecule  with  the  resultant  sharing  of  a 
hydrogen  proton.  Each  free  molecule,  therefore,  has  tiie  capability  of  having 
four  neighboring  molecules  bonded  to  it.  In  fact,  this  model  might  predict  that 
all  of  the  HjO  molecules  might  be  bonded  together.  Of  course  this  cannot  be 
so  because  of  the  liquid  nature  of  water.  This  is  not  even  true  for  ice  (Hasted 
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1973).  Further  experimental  evidence  is  needed  to  help  postulate  a  reasonable 
pure  liquid  water  physical  structure. 

Given  that  the  individual  water  molecules  are  dipoles  (which  is  supported 
by  the  fact  tliat  pure  water  does  exhibit  a  relatively  high  static  jjermittivity  of 
around  80),  then  one  could  reasonably  assume  that  the  dominant  high- 
frequency  loss  mechanism  showTi  in  Figure  1  is,  in  fact,  associated  with  the 
dipoles  not  being  able  to  keep  up  with  the  rapidly  changing  electric  field. 
Switching  off  a  static  electric  field  applied  to  a  sample  of  water  would  produce 
a  finite,  albeit  small,  delay  in  the  material  returning  to  an  unpolarized  state 
(relaxation).  Such  behavior  could  be  modeled  by 


P  a  e 


(11) 


where 

/*  =  the  polarization  of  the  material 
r  =  time 

X  =  a  characteristic  of  the  material  called  its  relaxation  time 
The  rate  of  change  of  polarization  is  then 


dP 

~di 


a 


1 
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(12) 


Carrying  this  argument  a  bit  further,  one  can  assume  that  the  relaxation 
process  behaves  as  a  temperature  dependent  chemical  reaction  for  which 


rale  of  reaction  -  k  s] 


(13) 


where 

/k  =  the  rate  constant,  or  the  probability  of  the  reaction  taking  place 
[S]  =  the  concentration  of  reacting  material 
Now,  Arrhenius  postulated  that  the  rate  constant  goes  like  (Grimshaw  1971) 
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(14) 


^  A 

k  ^  Ce 


where 

i4  =  the  activation  energy,  or  the  energy  needed  to  overcome  some 
equilibrium  state 

/?  =  the  universal  gas  constant 

r  =  the  absolute  temperature 

Drawing  parallels  between  the  hypothesized  time  rate  of  change  of  material 
polarization  and  Arrhenius'  rate  of  chemical  reaction  by  thinking  of  polariza¬ 
tion  as  concentration  and  the  inverse  of  the  relaxation  time  as  the  rate  constant, 
one  could  write  that 


J.  =  Ce  (15) 

t 


or 


In 


A 

fO 

a  -  _ 

R 

(16) 


If  measurements  of  y/t  (model  developments  in  the  next  section  will  show  that 
I  lx  is  the  radial  frequency  at  peak  loss)  as  a  function  of  temperature  plotted 
semil  igarithmically  as  a  straight  line,  then  one  could  estimate  the  activation 
energy  required  for  that  process  to  take  place  from  the  slope  of  the  curve. 

Ray’s  empirical  fit  for  liquid  water  includes  a  relationship  between  the 
relaxation  wavelength  and  temperature  that  looks  just  like  Equation  15.  A 
simple  calculation  results  in  an  estimate  for  activation  energy  of  water  to  be 
about  5.4  kcal/mole,  which  is  about  the  energy  required  to  break  a  hydrogen 
bond  (Cotton  and  Wilkinson  1972).  In  other  words,  hydrogen  bond  breaking 
becomes  a  likely  candidate  for  describing  the  relaxation  process  in  liquid  water 
and  is,  therefore,  an  indication  of  water  siruciure. 
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Debye  relaxation  model 

The  experimental  data  on  the  electrical  behavior  of  water  clearly  demon¬ 
strate  a  frequency-dependent  response  or  an  anomalous  dispersion  within  the 
1-  to  100-GHz  frequency  range.  The  preceding  section  indicated  that  the 
structure  of  liquid  water  includes  a  collection  of  permanent  dipoles;  i.c..  pani¬ 
cles  subject  to  realigitment  in  the  presence  of  an  oscillating  electromagnetic 
field.  An  obvious  analogy  found  in  classical  mechanics  for  an  oscillating 
particle  whose  behavior  changes  with  frequency  is  the  harmonic  oscillator. 
What  follows  is  the  development  of  a  simple  model  for  the  dipole  relaxation 
loss  seen  in  liquid  water  based  on  the  harmonic  oscillator. 

Using  the  notation  and  following  (with  some  extension)  the  treatment  of 
optical  dispersion  in  materials  presented  by  Reitz.  Milford,  and  Christy  (1980), 
one  can  consider  the  dipole  to  behave  like  a  one-dimensional,  damped,  forced 
oscillator  where  the  positive  charge  e  moves  and  the  negative  charge  remains 
fixed. 


m  -t-  G  +  Cx  =  eE  (17) 

dt^  dt 


where 

X  =  a  measure  of  the  displacement 
G  =  a  viscosity  coefficient 
C  =  a  spring  resistance  coefficient 
e  =  the  value  of  charge  at  each  end  of  the  dipole 
=  the  local  electric  field 

If  the  driving  force  were  zero  and  the  particle  were  given  some  initial  displace¬ 
ment,  x^,  then 


Cl  _t 

x(r)  =  XqC  ^  =  XqC  ^ 


(18) 


which  is  consistent  with  the  previous  description  of  relaxation  behavior  and 
activation  energies. 

Now  assume  that  both  the  driving  force  (the  local  electric  field)  and  the 
resulting  displacement  vary  sinusoidally  with  amplitudes  Ap  and  A^, 


Chapter  2  Electrical  Propertes  of  Water 


respectively,  and  that  the  inertia  forces  are  negligible.  One  can  then  write  the 
equation  of  motion  as 


G  (-iOVl.e  +  CA^e  =  eA^ 


(19) 


or 


= 

'  C  -  iGoi 


(20) 


For  the  single  dipole,  the  dipole  moment  is  eA^.  Therefore,  the  polarization  of 
the  material  P,  which  is  the  number  of  dipole  moments  per  unit  volume  and  is 
proportional  to  the  amplitude  of  the  driving  force  through  the  suscepiability  X- 
goes  like 


P  =  NeA^  = 


(21) 


where  =  the  number  of  dipt'les/unii  volume. 
Then  one  has 


X  = 


Ne^ 

~ 

1  -  1(01 


(22) 


which  shows  that  the  electric  susceptibility  of  the  medium  is  both  complex  and 
frequency  dependent.  But  to  model  this  behavior  in  terms  of  the  dielectric 
constant,  one  must  then  pose  the  relationship  between  dielectric  constant  and 
susceptibility  as 


e  =  e_  +  4rcx 


(23) 


where  £_  replaces  the  unity  term  in  the  usual  definition  for  linear  dielectrics 
and  represents  the  high-frequency  limit  on  the  real  pan  of  the  dielectric  con¬ 
stant.  Combining  the  last  two  expressions,  one  can  now  write 
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e  -  £ 


AriNe^ 

C 


(24) 


1  -  /on 


Referring  to  the  static  (to  =  O)  dielectric  constant  as  £„.  one  then  has  that 


AnNe^ 

“T~" 


=  Co  ~ 


(25) 


or  finally. 


e  -  e. 


=  ^  ~ 

1  -  /  on 


(26) 


Recalling  that 

ic" 


the  above  expression  can  be  solved  for  e'  and  e". 


e'  =  e_ 


1  +  (on)^ 


(27) 


1  +  (toxf 


(28) 


These  are  the  often-referenced  Debye  equations  (Debye  1929)  for  modeling  the 
dielectric  behavior  of  materials  made  up  of  polar  molecules.  Plotted  on 
semilogarithmic  scales,  these  equations  approximate  the  anomalous  behavior 
shown  for  liquid  water  on  Figure  1.  At  on  =  1,  e"  is  maximized,  which  says 
that  the  maximum  loss  due  to  the  dielectric  relaxation  mechanism  occurs  at  a 
frequency  equal  to  the  inverse  of  the  relaxation  lime. 

Eliminating  on  from  these  equations  allows  one  to  write  the  expression 
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+  :t"  j  ‘  = 

Co  -  E. 

(29) 


Thus,  the  Debye  equations,  '.vhen  drawn  in  e',  e"  space,  result  in  a  circle  cen¬ 
tered  at 


2 


e"  =  0 


and  having  a  radius  of  (e^  ‘  O/  2.  Such  a  plot,  shown  in  Figure  3,  is  known 
as  a  Cole-Cole  diagram.  Three  different  relaxation  frequencies  representing 
the  three  temperatures  shown  in  the  figure  were  chosen  using  Ray's  empirical 
model.  One  should  further  note  that  at  to  =  l/x  (the  inverse  of  the  relaxation 
time) 

2  2 


and 


f." 


2 


i.e..  the  peak  relaxation  loss  is  represented  by  the  point  at  the  top  of  the  Cole- 
Cole  diagram. 

What  has  been  demonstrated  in  the  preceding  paragraphs  is  u'lai  a  simple 
mechanical  model  can  do  a  very  good  job  in  representing  the  high-frequency 
dispersive  dielectric  behavior  of  liquid  water.  Of  course,  a  simple  electrical 
analog  can  also  be  developed.  In  terms  of  simple  circuit  elements,  the  equiva¬ 
lent  circuit  that  precisely  models  the  Debye  equations  is  shown  in  Figure  4 
(Cole  and  Cole  1941)  where  e  =  e'  -t-  re"  is  the  equivalent  capacitance  of  the 
circuit  (Equation  26).  Simple  models  like  this  will  prove  to  be  very  helpful  in 
analyzing  the  behavior  of  the  complex  dielectric  constant  in  moist  soils. 


Cole-Cole  relaxation  model 

Cole  and  Cole  (1941)  observed  that  a  good  bit  of  experimental  data  on  both 
polar  liquids  and  polar  solids  were  not  fit  by  the  semicircular  e'.  c" 
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Figure  3.  Graphical  representation  of  Debye  equations  in  e',  e'  space 


space  predictions  of  Debye’s  equations.  Rather,  it  seemed  that  over  large  fre¬ 
quency  intervals,  the  data  seemed  to  be  best  fit  by  circular  arcs;  i.e.,  pieces  of 
circies  in  e',  e"  space  that  were  centered  below  the  e"  =  O  line.  Using  a  series 
of  geometrical  and  analytical  arguments,  the  authors  showed  that  linear  materi¬ 
als  that  behaved  in  this  manner  can  be  described  by  the  relationship 


£  -  e. 


Eo  -  £- 

1  -  {imr 


(30) 


Expressions  for  e'  and  e"  found  by  substituting 
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Figure  4.  Equivalent  circuit  tor  Debye  rnoctel 


are  those  slightly  modined  by  Ray  (1972)  to  produce  the  curves  in  Figure  1. 

Furthenr.ore,  it  is  clear  from  the  form  of  Equation  30  that  a  simple  equiva¬ 
lent  circuit  can  still  be  drawn  for  materials  that  follow  the  behavior  described 
above.  Figure  5  shows  that  the  resistance  in  the  Debye  equivalent  circuit  now 
becomes  a  complex  impedance.  Representative  plots  of  the  Cole-Cole  relax¬ 
ation  model  arc  shown  on  Figure  6.  These  curves  were  derived  from  Ray's 
empirical  model  for  which  artificially  large  alpha  values  were  selected.  The 
venical  lines  on  the  figure  reflect  the  contribution  of  the  conductivity  term 
added  to  the  Cole-Cole  model  by  Ray  and  were  not  pan  of  the  original  devel¬ 
opment.  They  do,  however,  give  an  indication  of  the  way  real  material 
responses  may  look. 
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Figure  6.  Representative  Cole-Cole  plots 
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3  Electrical  Properties  of 
Moist  Soils 


A  quick  review  of  available  dieleciric  data  on  dry  soils  and  wet  soils  leads 
to  the  conclusion  that  the  dominant  factor  that  controls  the  electrical  behavior 
of  soils  is  the  presence  of  water.  Over  a  broad  frequency  range,  the  real  pan 
of  the  complex  relative  dielectric  constant  (referred  to  in  this  text  as  permittiv¬ 
ity)  of  dry  soil  minerals  changes  very  little,  covering  a  range  of  values  from 
about  2  to  about  6  (Nelson,  Lindroth,  and  Blake  1989;  Ulaby  et  al.  1990). 
However,  when  considered  from  the  perspective  of  the  idealized  elearomag- 
netic  wave  reflection  phenomenon  for  lossless  materials  as  described  in  Appen¬ 
dix  B,  this  can  translate  into  reflection  coefficients  that  vary  from  -0.17  to 
-0.42.  In  terms  of  power  (which  is  proportional  to  radar  backscatter  coeffi¬ 
cients).  these  small  variations  in  dry  soil  properties  result  in  reflectances  that 
span  values  from  0.029  to  0.176,  or  a  range  of  about  8  db. 

When  water  is  added  to  the  soil  fabric,  further  substantial  changes  in 
reflectance  at  the  soil/air  interface  can  take  place  Take  the  permittivity  of 
water  to  be  about  80  (at  low  frequencies).  Then  for  normal  incidence  electro¬ 
magnetic  waves  and  lossless  media,  one  can  calculate  a  reflection  coefficient 
of  about  -0.8  and  a  reflectance  of  about  0.64.  Now,  take  the  permittivity  of  a 
very  moist  soil  to  be  about  20.  Tlien  under  ideal  conditions,  its  reflectance  at 
normal  incidence  could  be  as  much  as  0.4.  In  other  words,  a  very  moist  soil 
could  theoretically  reflect  more  than  thineen  times  the  power  (or  more  than 
1 1  db)  than  can  the  driest  soils.  While  airborne  scatterometers  may  easily  see 
this  much  variation  in  measurements  over  large  areas  of  so-called  homo¬ 
geneous  terrain  due  to  surface  roughness  effects  and  interference  phenomena, 
these  differences  in  electrical  properties  of  soils  could  be  quite  significant  for 
close-up  measurements  with  ground-based  electromagnetic  devices. 

As  an  example  of  the  effect  of  soil  moisture  on  real-world  situations,  con¬ 
sider  the  data  reported  by  Ulaby,  Cihlar,  and  Moore  (1974)  that  was  collected 
with  a  ground-based  radar  scatterometer  located  in  an  unplantcd,  plowed  field 
of  clay  loam  soil.  As  shown  in  Figure  7,  the  power  returned  to  the  scatiero- 
meter  while  looking  nearly  straight  down  at  the  soil  surface  (0-deg  incidence) 
increased  by  20  db  as  moisture  contents  increased  from  4.3  to  36.3  percent 
(average  values  in  the  first  5  cm).  Twenty  decibels  means  a  lOO-fold  increase 
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Figure  7.  Radar  backscatter  coefficient  versus  moisture  content  (from  Ulaby, 
Cihlar,  and  Moore  (1974)) 


in  power.  The  arguments  above  would  indicate  an  anticipated  bound  of  about 
11  or  12  db  increase.  It  is  possible  that  surface  roughness  effects  and/or  con¬ 
structive  interferences  could  account  for  the  larger  field  measurement  numbers. 
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Data  on  Variations  in  Frequency,  Moisture,  and 
Temperature 


At  this  point,  it  would  be  pnident  to  assess  what  is  already  known  about  the 
electrical  behavior  of  moist  soils.  What  are  the  conditions  under  which  data 
have  already  been  collected?  What  kinds  of  soils  have  been  studied?  Is  there 
general  agreement  on  results  of  previous  measurements?  Have  the  measure¬ 
ments  led  to  any  insight  as  to  the  loss  mechanisms  involved? 

While  the  governing  equations  for  the  propagation  of  electromagnetic 
waves  in  linear,  homogeneous,  isotropic  media  have  been  accepted  since  the 
late  1800s,  only  in  the  last  30  years  or  so  has  any  serious  attention  been  given 
to  the  complex  dielectric  response  of  heterogeneous  mixtures  such  as  soils. 
Nevertheless,  all  measurements  of  soil  electrical  properties  are  still  interpreted 
in  a  macroscopic  sense  as  if  the  soil  is  truly  homogeneous  and  isotropic. 

Table  1  identifies  some  of  the  most  relevant  work  on  measuring  the  electrical 
response  of  materials  that  ultimately  apply  to  this  study  on  the  response  of 
soils.  No  such  list  could  ever  hope  to  be  complete,  but  it  does  serve  to  pro¬ 
vide  a  historical  perspective  and  a  proper  point  of  departure  for  the  work  to  be 
conducted  in  this  and  future  studies.  Paragraphs  that  follow  will  highlight 
several  of  these  contributions. 


Structure  of  clay  minerals 

Although  this  section  may  seem  out  of  place,  it  is  important  to  set  the  stage 
for  discussions  of  soil  electrical  behavior  by  briefly  addressing  the  structure  of 
the  soil  elements  that  probably  are  most  influential  in  determining  their  electri¬ 
cal  response;  namely,  the  clay  minerals. 

Day  minerals  are  layered  silicates  whose  fundamental  building  blocks  are 
tetrahedral  sheets  in  which  the  teirahedra  are  linked  at  their  comers  and  octa¬ 
hedral  sheets  in  which  the  ociahedra  are  linked  along  their  edges  (Moore  and 
Reynolds  1989).  Each  tetrahedron  is  formed  from  four  oxygen  ions  normally 
surrounding  Si‘*  cations  (which  can  also  be  replaced  by  Al’*  or  Fe^'"  cations). 
Each  octahedron  is  formed  from  six  oxygen  ions  (or  hydroxyl  ions)  surround¬ 
ing  a  cation  that  is  normally  either  Al’*,  Mg^*,  Fe^*,  or  Fe’*. 

Most  clay  minerals  fit  into  structural  classifications  referred  to  as  1:1  or 
2:1.  A  1 :1  structure  is  comprised  of  a  teuahedral  sheet  joined  to  an  octahedral 
sheet.  The  mechanism  for  this  bonding  is  that  the  apical  oxygens  of  the  tetra- 
hedra  replace  two  out  of  every  three  anions  in  the  octahedral  sheet.  The 
2:1  stmeture  is  generated  by  a  second  tetrahedral  sheet  bonding  to  the  opposite 
side  of  the  octahedral  sheet.  Slight  differences  in  sheet  dimensions  (or  anion 
.spacing)  due  to  various  combinations  of  cations  can  cause  distortions  of  the 
layered  structure,  sometimes  so  severe  as  to  result  in  tubular  geometries. 
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Clay  mineral  layers  formed  by  the  tetrahedral  and  octahedral  sheets  are 
sometimes  electncally  neutral,  but  most  often  are  somewhat  negatively  charged 
due  to  the  substitution  of  lesser  valence  cations  for  Si^'  in  the  tetrahedral  sheet 
and/or  Al’*  in  the  octahedral  sheet.  Furthermore,  because  of  the  finite  lateral 
dimensions  of  clay  mineral  crystals,  unsatisfied  bonds  exist  at  the  edges  of  the 
layers  tfiat  also  result  in  layer  charge  imbalance. 

The  following  simplified  descriptions  of  the  clay  minerals  that  will  be 
referred  to  in  this  study  can  now  be  given.  Kaolinite  is  a  1:1  structure  mineral 
which,  if  in  a  very  pure  form,  will  have  little  or  no  layer  charge  due  to  cation 
exchange,  but  will  attract  ions  or  polar  molecules  such  as  water  to  its  edges. 
lUitc,  montmorillonite,  and  hectorite  arc  2;1  minerals.  Many  opinions  seem  to 
exist  regarding  what  illite  really  is  in  terms  of  a  structural  formula,  but  it  is 
generally  accepted  that  it  is  a  mica  material  whose  layers  have  a  half-unit-cell 
charge  imbalance  of  about  0.8  and  whose  interlayer  spaces  are  occupied  by 
cations  (usually  potassium)  whose  spatial  distribution  nearly  balances  the  layer 
charges.  Montmorillonite  and  hectorite  have  smaller  layer  charge  imbalances; 
the  interlayer  cations  can  attract  water,  which  results  in  a  swelling  of  the  lay¬ 
ered  structure  and  the  creation  of  hydrogen  bonds  between  the  water  molecules 
and  the  anions  of  the  tetrahedral  sheet  surfaces  (Moore  and  Reynolds  1989). 
These  swelling  clays  have  a  greater  affinity  for  water,  which  should  cause 
much  different  electrical  responses  than  for  the  nonswelling  kaolinite. 


Dispersive  behavior  in  moist  soiis 

Smith  (1971)  conducted  a  series  of  tests  on  the  electrical  behavior  of 
saturated  clays  up  to  a  maximum  frequency  of  about  60  MHz  using  a  capaci¬ 
tive  bridge  measurement  apparatus.  He  looked  at  three  distinct  clay  types:  a 
montmorillonite,  an  illite,  and  a  kaolinite. 

A  representative  set  of  Smith’s  measurement  results  is  given  in  Figure  8  in 
terms  of  the  real  part  of  the  dielectric  constant  and  conductivity  versus  fre¬ 
quency.  Of  particular  interest  is  the  observation  that  clays  with  higher  basic 
structural  water  content  have  larger  complex  dielectric  constants  than  those 
with  less  water.  Conclusions  should  not  be  drawn  from  consolidation  differ¬ 
ences.  Density  should  clearly  affect  the  dielectric  properties  of  the  basic  min¬ 
eral  themselves  by  providing  varying  degrees  of  panicle  contact  and,  hence, 
varying  electrical  path  lengths  through  the  matrix;  however,  the  degree  of 
effect  is  not  currently  understood. 

Lundien  (1971)  reported  permittivity  and  conductivity  measurements  over  a 
range  of  frequencies  between  10  MHz  and  1.5  GHz.  At  the  lowest  frequen¬ 
cies,  he  used  a  capacitive  bridge  setup.  Data  at  about  300  MHz  were  collected 
from  radar  reflectance  measurements.  The  higher  frequency  data  were 
collected  using  what  was  referred  to  as  a  microwave  interferometer,  but  was, 
in  fact,  a  free-space  transmission  measurement  apparatus  with  the  sample  faces 
tilted  to  the  radar  path. 
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Figure  8  Dielectric  constant  ol  clays  versus  frequency  (from  Smith  (1971))  (Sheet  1  ol  3) 
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Figure  8.  (Sheet  3  of  3) 


Lundien’s  results  are  shown  in  Figure  9  for  three  types  of  soil;  a  very 
poorly  graded  sand,  a  well-graded  silt,  and  a  high  clay  content  locally  available 
soil  (50  percent  of  particles  by  weight  <.0075  mm).  Conclusions  are  difficult 
to  draw,  because  the  data  include  a  broad  range  of  moisture  levels.  The  Long 
Lake  clay  is  described  by  Lundien  as  being  composed  mostly  of  montmorillon- 
ite  panicles. 

A  ver>'  nice  set  of  moist  soil  electrical  measurements  over  a  broad  range  of 
frequencies  was  collected  by  Hockstra  and  Delaney  (1974)  using  a  7-mm  coax¬ 
ial  line  for  frequencies  up  to  about  3  GHz  and  a  series  of  slotted  waveguide 
devices  from  about  5  GHz  to  about  20  GHz.  Copies  of  their  reponed  data  as  a 
function  of  frequency  are  shown  in  Figures  10  and  11  for  two  different  clays. 
Figure  10  (a)  shows  data  for  one  gravimetric  moisture  content  and  two  temper¬ 
atures.  while  Figure  10  (b)  shows  data  for  temperature  and  two  moisture  con¬ 
ditions.  Qualitatively  they  observed  a  Debye-like  dispersion  in  the  moist  soils, 
but  at  a  lower  frequency  than  for  bulk  liquid  water.  They  also  observed  that 
increasing  temperatures  reduced  the  frequency  of  maximum  dielectric  loss 
(which  is  inconsistent  with  the  response  of  water)  as  did  increasing  water 
content.  They  further  attempted  to  fit  a  modified  Debye  relaxation  equation  to 
the  data  fo  one  of  the  clays  (Figure  11)  with  partial  success.  They  concluded 
that  there  was  no  difference  between  the  relaxation  of  water  in  sandy  soil  and 
that  for  clays,  which  contradicts  the  thesis  that  chemical  bonding  of  the  water 
dipoles  to  clay  particles  should  cause  a  shift  of  the  relaxation  frequency. 

Hallikaincn  et  J.  (1985)  collected  dispersion  data  using  waveguide  trans¬ 
mission  techniques  in  the  1-  to  2-  and  4-  to  6-GHz  bands  and  a  free  space 
transmission  technique  at  eight  selected  frequencies  between  4  and  18  GHz. 
They  conducted  measurements  on  soil  types  ranging  from  sand  loam  (51  per¬ 
cent  sand  and  13  percent  clay)  to  a  silty  clay  (5  percent  sand  and  47  percent 
clay).  Figure  12  contains  their  published  results  for  one  soil  (42  percent  sand. 
49.5  percent  silt,  and  8.5  percent  clay)  collected  at  an  ambient  temperature  of 
about  23  ”C. 

Anouher  set  of  fretiuency-dependent  data  is  that  collected  by  Campbell 
(1988)  using  a  tined  coaxial  probe  and  a  device  for  precisely  controlling  sam¬ 
ple  moisture  content  and  sample  temperature.  Frequency  was  limited  to 
50  MHz  for  these  measurements.  The  objective  of  these  studies  was  to  exam¬ 
ine  the  low-frcquency  dielectric  loss  mechanisms  in  moist  soils. 

Figure  13  is  a  summary  of  Campbell's  measurements  for  one  type  of  silty 
soil  that  shows  the  low-frequency  dispersion  of  this  soil  as  a  function  of  volu¬ 
metric  water  content.  Obviously,  one  is  not  able  to  ascenain  the  changes  in 
Debye  tvpc  relaxation  because  of  the  frequency  limitations  of  these  data.  Of 
panicular  note  was  the  conclusion  drawn  by  Campbell  that  the  predominant 
loss  mechanism  in  this  frequency  range  was  that  of  ionic  conductivity  and  not 
the  Maxwell- Wagner  effect  h>pothcsizcd  by  Smith  (1971).  These  mechanisms 
Will  be  discussed  in  a  later  section. 
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Dielectric  constant  of  loam  versus  frequency  (from  Hallikainen  et  al.  (1985)) 


IMAGINARY  WEtECTRIC  CONST. 


Figure  13.  Dieleciric  constant  of  Manchester  silt  (from  Campbell  (1988)) 
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Existing  data  on  effects  of  moisture  content 


Because  of  the  obvious  applications  of  remote  sensing  technology  to  the 
determination  of  soil  moisture  levels,  considerable  data  have  been  collected 
and  reported  in  the  literature.  An  early  example  is  that  of  Lundien  (1966)  in 
which  he  used  four  different  radar  scattcrometers  and  large  specially  prepared 
soil  samples  to  measure  the  apparent  dielectric  constant.  “Apparent”  ip  this 
context  means  that  the  dielectric  constant  was  computed  from  the  square  of  the 
Fresnel  coefficients  (Appendix  B)  assuming  no  imaginary  components  (no 
losses).  Lundien’s  results  for  two  soil  types  are  shown  in  Figures  14  and  15  in 
which  he  suggests  an  exponential  fit  might  be  in  order.  The  data  are  limited 
in  quantity  because  of  the  difficulty  of  sample  preparation  and,  as  such,  do  not 
reveal  what  later  researchers  believe  to  be  a  bilinear  type  of  response. 

Lundien’s  later  work  (1971)  using  free  space  transmission  measurement 
techniques  is  summarized  in  Figures  16.  17,  and  18.  The  moisture  content 
reported  here  is  volumetric  (volume  of  waterAolume  of  sample)  as  opposed  to 
his  earlier  gravimetric  (weight  of  water  in  sample/dry  weight  of  soil  in  sample; 
rti,  =  data.  Polarization  labels  on  these  figures  refer  to  the  orien¬ 

tation  of  the  antennas  with  regard  to  the  horizon.  The  data  show  no  polari¬ 
zation  effects,  even  for  the  clay,  which  under  compaction  can  assume  a 
somewhat  regular  platelet  stmeture  or  fabric.  The  quantity  of  data  collected  is 
much  more  substantial  than  the  previous  results  reported  above,  and  in  some 
cases,  shows  a  bilinear  response.  One  can  argue  from  these  data  that  soil 
fabric  does  not  play  a  major  role  in  detennining  the  complex  dielectric 
response  of  moist  soils. 

Hoekstra  and  Delaney  (1974)  also  found  little  difference  in  the  dielectric 
response  of  their  sand,  silt,  and  clay  samples  as  shown  in  the  data  on  Fig¬ 
ure  19.  From  an  engineering  perspective,  their  conclusion  is  valid.  However, 
if  one  wishes  to  examine  the  physical  mechanisms  behind  the  electrical  losses 
in  moist  soils,  one  would  have  to  have  more  data  on  each  soil  type  than  is 
shown  in  Figure  19. 

Hallikainen  et  al.  (1985)  reported  results  of  measurements  on  five  different 
soil  types.  Figure  20  shows  curves  of  dielectric  constant  versus  volumetric 
water  content  drawn  from  polynomial  fits  to  the  available  data.  One  could 
easily  argue  that  these  results  indicate  some  differences  in  electrical  behavior 
due  to  soil  texture  at  relatively  low  frequencies  (<5  GHz),  but  hardly  any  dif¬ 
ferences  at  higher  frequencies.  Their  assumption  of  a  polynomial  fit  (second 
order  in  volumeuic  moisture  content)  precludes  any  discussion  of  bilinear 
behavior. 

Another  source  of  experimental  data  on  moisture  content  variation  is  that 
reported  by  Campbell  (1988)  and  shown  in  Figure  21  for  three  different  soils. 
Clearly,  at  low  frequencies  (<  50  MHz)  dielectric  energy  loss  goes  up  with 
moisture  content  and  is  very  dependent  upon  soil  type  at  the  upper  end  of  the 
frequency  range  of  these  experiments. 
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Figure  14.  Apparent  dielectric  constant  of  a  silt  loam  versus  gravimetric  moisture  (from  Lundien 
(1966)) 


Effects  of  sample  temperature  on  dielectric  properties 

While  several  of  the  preceding  figures  contain  information  on  the  effects  of 
sample  temperature  on  dielectric  propenies.  special  attention  should  be  given 
to  the  work  done  by  Hoekstra  and  Delaney  (1974).  While  they  were  not  able 
to  conclude  anything  on  the  temperature  dependence  of  the  characteristic  fre¬ 
quency  of  Debye-type  relaxation,  they  did  rcpon  on  some  very  interesting 
dielectric  behavior  near  the  freezing  point  of  water.  Figure  22  shows  how 
frozen  soils,  in  general,  have  a  lower  dielectric  response  than  do  unfrozen 
soils,  and  that  increasing  moisture  content  reveals  some  very  anomalous 
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Figure  15  Apparent  dielectric  constant  of  a  clay  versus  gravimetric  moisture  (Irom  Lundien 

(1966)) 

behavior  near  0  “C.  Some  earlier  results  reported  by  Hocksira  and  Doyle 
(197 1)  support  the  same  contcr.iion. 

Another  very  informative  study  was  carried  out  by  Delaney  and  Arcone 
(1982)  that  involved  the  cross-plotting  of  extensive  lime  domain  rctleclomeiry 
data  to  produce  the  curves  shown  in  Figure  23  for  a  sill  and  a  sand. 

Hallikaincn  ct  al.  (1985)  also  reported  measuremenus  made  at  four  different 
sample  temperatures  ranging  Irx^m  -^23  to  -24  °C.  Their  resulus,  showTi  in 
Figure  24,  reveal  the  same  relative  insensitivity  to  temperature  at  subfrec/ing 
conditions  as  was  reported  by  others. 
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Figure  17.  Dieleciric  constant  ct  a  silt  versus  volumetric  moisture  (from 
Lundien  (I97i)) 
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Figure  18.  Dielectric  constant  of  a  clay  versus  volumetric  moisture  content 
(from  Lundien  (1971)) 
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Figure  19.  Dielearic  constant  of  soils  versus  volumetric  moisture  at  10  ”0  (from  Hoekstra  and 
Delaney  (1974)) 


Bound  versus  Free  Water 

In  spite  of  conflicting  conclusions  regarding  the  dependence  of  soil  dielec¬ 
tric  response  on  soil  texture,  many  researchers  contend  that  the  dependetKc  is 
there  CWang  and  Smugge  1980). 

If  moist  soils  with  different  textures  do,  in  faa.  possess  different  dielectric 
properties  for  the  same  levels  of  moisture  contenu  it  must  be  because  of  differ¬ 
ent  soil-water  bonding  mechartisms.  Qay  particles  offer  a  relatively  large 
amount  of  electrostatically  charged  surface  area  (van  Olphen  1963;  IXrbson  ct 
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Figure  20.  Dielectric  constant  ot  soils  versus  volumetric  rTX)isture  (from  Hallikainen  et  al. 
(1985)) 
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Figure  21 .  Dielectric  constant  tor  soils  versus  volumetric  moisture  (trorr. 
Campbell  (1988))  (Sheet  i  of  3) 


al.  1985)  that  could  support  numerous  bonding  mechanisms.  The  unique  smic- 
ture  of  clayey  soils  also  present  many  opponuniiies  to  alter  the  dielectric 
response  of  the  soil  as  a  whole  by  the  creation  of  microscopic  capacitive  ele¬ 
ments.  Even  if  the  surface  charge  density  of  sand  particles  is  not  significantly 
different  from  that  of  clay  particles  (van  Olphen  1963),  the  rounded  shape  of 
the  sand  particles  and  their  smaller  specific  surface  area  (square  meters/gram) 
must  lead  to  different  mechanical  as  well  as  electrical  behavior. 
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Figure  21  (Sheet  2  of  3) 

Relating  the  physical  behavior  of  moist  soil  to  its  dielectric  response  may 
even  be  possible.  For  example,  soil  scientists  have  long  associated  varying 
levels  of  soil  tension,  reflected  in  the  pressure  required  to  force  water  from  a 
sample,  with  varying  degrees  of  soil  water  bonding.  Wang  and  Smugge 
(1980)  took  the  concept  of  soil  tension  a  bit  further  by  correlating  a  particular 
level  of  tension,  known  as  the  wilting  point,  to  something  called  the  transition 
moisture  content  in  moist  soils.  The  uansition  moisture  content  is  the  break 
point  in  the  bilinear  dielectric  response  referred  to  in  an  earlier  section. 
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Universal  agreement  on  a  single  mechanism  for  soil-water  bonding  does  not 
exist.  Rather,  several  hypotheses  are  offered,  each  of  which  can  be  both  sup¬ 
ported  aiid  refuted  by  experimental  evidence  (Mitchell  1974). 
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(b)  North  Slope  sand 


Figure  23.  Dielectric  constant  of  soils  versus  temperature  at  several  volumetric 
moisture  contents  and  a  frequency  of  0.5  GHz  (from  Delaney  and 
Arcone  (1982)) 
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Dielectric  constant  for  two  soils  as  a  function  of  frequency  and  temperature  (from  Hallikainen  et  al.  (1985)) 


Hydrogen  bonding 


The  most  easily  supportable  picture  of  soil-water  bonding  is  that  of  hydro¬ 
gen  bonding,  which,  as  discussed  previously,  is  the  sharing  of  a  proton  by  two 
electronegative  atoms.  Mitchell  (1976)  points  out  that  soil  particle  surfaces  are 
usually  composed  of  either  a  layer  of  oxygen  atoms  or  a  layer  of  hydroxyl 
units.  He  argues  that  the  shared  proton  is  supplied  by  the  water  molecule  if 
the  surface  layer  is  oxygen  atoms,  while  the  hydroxyl  supplies  the  proton  in 
the  other  situation.  In  either  case,  activation  energy  measurements  do  provide 
some  supporting  evidence  of  this  mechanism.  Newman  (1987)  further  explains 
that  unsatisfied  bonds  at  the  edges  of  clay  panicles  present  further  bonding 
opportunities  to  the  water  molecules  depending  on  the  pH  of  the  fluids  in  the 
moist  soil. 


van  der  Waal's  forces 

A  second  proposed  mechanism  for  soil-water  bonding  is  that  of  electrostatic 
attraction  combined  with  the  dipole  nature  of  the  water  molecule  to  cause 
several  layers  of  molecules  to  migrate  to  the  negatively  charged  surfaces  of  the 
soil  panicles.  In  clayey  soils  with  platelet  surfaces  nearly  parallel,  creating 
microscopic  dielectric-filled  capacitive  elements  through  the  van  der  Waal 
attraction  mechanism  is  possible  if  one  also  has  cations  available  to  complete 
the  dipole  path  as  shown  in  Figure  25. 


Hydration  of  exchangeable  cations 

A  third  soil-water  bonding  mechanism  is  a  parasitic  one  in  which  cations 
are  attracted  to  the  negatively  charged  particle  surfaces.  Newman  (1987) 
views  this  as  the  primary  water  adsorption  mechanism  in  swelling  clays.  This 
mechanism  assumes  that  the  cations  take  their  water  of  hydration  along  with 
them  for  the  ride.  Naturally,  one  must  have  cations  available  in  solution  for 
this  model  to  have  any  validity,  and  the  concentration  of  cations  determines 
the  amount  of  water  that  can  be  bonded.  Dobson  et  al.  (1985)  was  motivated 
by  such  a  mechanism  to  model  the  contribution  of  bound  water  in  a  four- 
component  dielectric  mixing  model.  In  reality,  a  double-layer  model  was  used 
for  ion  distribution  near  clay  particle  surfaces  to  moaify  bulk  conductivity  for 
free  water  losses  and  simply  assumed  values  of  the  complex  dielectric  constant 
for  bound  water. 

Osmosis 

The  fourth  mechanism  (and  that  least  supported  by  data)  is  that  due  to 
osmotic  forces.  Again  assuming  that  there  are  a  large  number  of  cations  in 
solution  and  that  their  attraction  to  the  negatively  charged  particle  surfaces  is 
stronger  than  the  attraction  of  the  dipole  water  molecules  to  the  surfaces,  this 
model  predicts  that  there  will  be  a  higher  concentration  of  cations  closer  to  the 
particle  surfaces.  The  existence  of  another  concentration  of  cations  on  a 
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Figure  25  Microscopic  capacitive  elements  in  clayey  soils 


nearby  surface  then  results  in  a  volume  of  lower  particle  concentration  between 
the  surfaces  and  the  resulting  migration  of  water  molecules  to  that  volume. 

Whatever  the  mechanism  is.  water  is  attracted  to  soil  panicles,  particularly 
clay  particles,  as  proven  by  the  extraordinary'  pressures  required  to  drain  soil 
samples.  In  the  iiumediair  vicinity  of  the  soil  panicle  surface,  the  concentra¬ 
tion  of  water  molecules  should  be  higher  than  in  the  bulk  or  free  water  occu¬ 
pying  the  void  spaces  (Martin  1960).  Thus,  one  has  to  believe,  even  if  the 
bonding  mechanism  was  that  of  hydrogen  bonds,  that  bound  water  should 
exhibit  a  different  dielectric  response  than  free  water  in  the  soil  matrix. 


Radio  Frequency  Loss  Mechanisms 

Clearly,  water  in  contact  with  soil  panicles  must  behave  differently  (electri¬ 
cally)  from  free  water.  This  shows  up  very  dramatically  at  radio  frequencies 
through  a  variety  of  loss  mechanisms  that  do  not  exist  for  the  dry  minerals  rr 
for  pure,  liquid  water.  Not  only  does  the  water  bind  to  the  soil  particle 
surfaces,  but  also  salts  are  dissolved  to  produce  a  source  of  conducting  ions 
(Mitchell  1974).  As  already  mentioned,  microscopic  capacitive  elements  can 
be  created  simply  from  the  plaielet-likc  structure  of  clay  panicles.  These  and 
other  phenomena  lead  to  a  number  of  loss  mechanisms  m  moist  soils  within 
the  radio  frequency  range  that  arc  discussed  in  the  following  paragraphs. 
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C  ,  ionic  conductivity 
DL  ,  charged  double  layers 
X  ,  crystal  water  relaxation 
I  ,  ice  relaxation 
MW  ,  Msxwell'Wagncr  effect 
S  surface  conductivity 
B  ,  bound  water  relaxation 
W1  ,  free  water  relaxation 


Figure  26.  Dielectric  loss  mechanisms  for  heterogeneous  moist  materials  (from  Hasted  (1973)) 

Hasted  (1973)  presented  a  particularly  useful  schematic  (Figure  26)  for  the 
relative  magnitude  and  frequency  range  of  application  of  each  of  what  are 
believed  to  be  the  major  contributors  to  dielectric  loss  in  heterogeneous  moist 
materials.  The  absolute  magnitudes  of  these  effects  should  not  be  regarded  as 
applying  to  any  particular  material.  A  very  brief  description  of  several  of 
these  mechanisms  follows.  More  detailed  expositions  may  be  found  in  other 
sources  (Campbell  1988;  Hasted  1973;  Mitchell  1976). 


Free  water  relaxation 

The  inability  of  free  water  dipolar  molecules  to  keep  up  with  an  applied 
alternating  field  has  already  been  discussed  in  Chapter  2.  Dipole  relaxation  in 
free  water  is  the  dominant  loss  mechanism  at  microwave  frequencies. 


Bound  water  relaxation 

Experimental  evidence  exists  (Muir  1954;  Hoekstra  and  Delaney  1974)  that 
indicates  a  lowering  of  the  critical  dipole  relaxation  frequency  in  water  that  is 
bound  to  solid  particle  surfaces.  One  can  think  of  this  physically  in  terms  of 
the  simple  mechanical  analog  described  in  Chapter  2  and  imagine  that,  for  all 
other  factors  kept  constant,  an  increase  in  the  viscosity  coefficient  will  result  in 
an  increase  in  the  relaxation  time  or  a  lowering  of  the  relaxation  frequency. 

Hall  and  Rose  (1978)  proposed  bound  water  relaxation  as  the  dominant  loss 
mechanism  in  carefully  prepared  kaolinite  clays  in  the  10-  to  l(X)-kHz  range  of 
the  spectrum,  but  their  conclusions  seem  to  contradict  the  observations  of 
others  and  the  implications  of  Figure  26. 
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Maxwell-Wagner  effect 


Another  loss  mechanism  often  referred  to  in  the  literature  as  the  dominant 
mechanism  at  radio  frequencies  is  one  termed  the  Maxwell-Wagner  effect 
(Campbell  1988;  Bidadi,  Schroeder,  and  Pinnavaia  1988).  It  originates  in  the 
experimental  observation  that  there  is  dispersion  in  a  suspension  of  conducting 
particles  in  nonconducting  dielectric  media  much  like  that  observed  in  a 
parallel -plate  capacitor  filled  with  two  different  dielectrics.  The  Maxwell- 
Wagner  effect  is  seen  as  an  accumulation  of  charge  at  the  interface  of  dissimi¬ 
lar  materials  during  the  flow  of  current  through  heterogeneous  material  due  to 
a  discontinuity  in  dielectric  constant  values.  This  charge  buildup  is  time 
dependent  and  results  in  dispersion  as  the  frequency  of  the  alternating  current 
increases  to  a  point  where  the  buildup  and  relaxation  cannot  keep  up. 

When  modeled  as  series  capacitors,  one  having  zero  conductivity  and  the 
other  a  nonzero  conductivity  leading  to  a  dielecuic  loss  through  Equation  10 
(assuming  the  pwlarization  loss  is  z-ro).  the  Maxwell-Wagner  effect  can  be 
described  by  Debye-like  relationships  (Hasted  1973:  Campbell  1988).  For 
example,  for  equal  thickr.ess  capacitors,  the  relaxation  time  expression  is  found 
to  be  as  follows; 


X 


47:Cj 


(31) 


Thus  the  relaxation  frequency  is  seen  to  get  smaller  as  conductivity  gets 
smaller. 

One  application  of  the  Maxwell-Wagner  model  (Bidadi,  Schroeder,  and 
Pinnavaia  1988)  had  water  as  the  conductive  medium  and  the  clay  platelets  as 
the  nonconductive  medium.  The  high  capacitance  of  the  closely  separated 
platelets  was  used  to  argue  for  the  observed  enhancement  of  clay-water  permit¬ 
tivity  above  the  level  for  pure  water. 


Surface  conductivity 

Another  radio  frequency  loss  mechanism  in  moist  heterogeneous  materials 
is  that  of  surface  conductivity.  A  pxjssible  scenario  for  this  loss  is  a  rapid 
intennolecular  proton  exchange  within  the  bound  water  layer  on  soil  particles 
(Hoekstra  and  Doyle  1971;  Fripiai  et  al.  1965;  von  Hippel  1988).  One  appli¬ 
cation  of  this  model  (Schwan  et  al.  1962)  led  to  the  calculation  of  a  Debye- 
like  relaxation  time 
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where  and  o,  refer  lo  the  permittiviiy  and  conductivity  of  the  aqueous  solu¬ 
tion  in  which  nonconducting  panicles  were  immersed.  Again,  as  with  the 
Maxweh-Wagner  effect,  reducing  the  conductivity  of  the  aqueous  solution 
reduces  the  relaxation  frequency. 

Charged  double  layers 

If  one  accepts  that  cations  can  be  attracted  to  negatively  charged  soil  parti¬ 
cle  surfaces,  then  one  can  po.ssibly  think  of  the  cations  surrounding  the 
charged  panicle  like  an  electron  cloud  surrounding  a  nucleus.  The  charged 
double  layer  loss  mechanism  then  has  a  parallel  in  the  concept  of  atomic  polar¬ 
izability.  As  the  electric  field  oscillates,  so  too  docs  the  ion  cloud  about  the 
panicle.  As  frequency  goes  up.  the  ability  of  the  cloud  to  diston  goes  down 
with  the  resulting  lowering  of  the  dielectric  loss 


Ionic  conductivity 

The  most  dominant  loss  mechanism  in  moist  heterogeneous  materials  at  low 
frequencies  is  that  due  to  ionic  conductivity,  the  4jio/a>  term  in  Equation  (10). 
Of  course,  this  conductivity  is  brought  about  by  dissolving  salts  in  the  soil  and 
providing  a  path  through  the  pore  spaces  for  the  transpon  of  electrical  energy 
to  take  place.  As  a  result,  low-frequency  ionic  conduebvity  in  moist  soils 
would  have  to  be  related  in  some  way  to  the  degree  of  saturation  or  how  the 
pore  spaces  fill  with  water. 


Activation  Energy  Data 

In  the  earlier  description  of  the  structure  of  water  found  in  Chapter  2  and 
under  the  section  entitled  Microwave  Frequency  Loss  Mechanisms,  arguments 
were  made  relating  the  frequency  of  peak  loss  in  a  Debye-type  dielectric  relax¬ 
ation  process  to  the  activation  energy  associated  with  that  process.  The  slope 
of  a  semilogarithmic  plot  of  that  relaxation  frequency  versus  the  inverse  of  the 
absolute  temperature  is  proportional  to  the  activation  energy. 

Hoekstra  and  Doyle  (1971)  argued  that  linear  plots  of  the  log  of  the  dielec¬ 
tric  loss  term  versus  the  inverse  of  absolute  temperature  can  also  be  obtained 
when  the  variation  in  dielccuic  loss  with  temperature  is  caused  by  changes  in 
the  thermal  energy  of  a  charge  carrier.  In  fact,  the  dynamical  theory  of  sorp¬ 
tion  (Hasted  1973)  relates  the  time  of  residence  of  a  water  molecule  on  a 
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surface  to  the  heat  of  adsorption  by  the  same  kind  of  exponential  relationship 
as  used  for  the  earlier  activation  energy  development.  Certainly,  if  loss  is 
inversely  proportional  to  this  residence  time  by  virtue  of  the  molecules  being 
free  to  reorient  in  the  alternating  electric  field,  then  Hoekstra’s  and  Doyle’s 
arguments  are  valid.  Thus,  if  the  dielectric  loss  at  any  given  frequency 
increases  with  sample  temperature,  the  semilogariihmic  plot  mentioned  above 
can  lead  to  an  estimate  of  the  energy  required  to  bring  about  such  a  change. 
This  change,  in  turn,  may  help  identify  the  loss  mechansims. 

In  their  1971  paper,  Hoekstra  and  Doyle  did.  in  fact,  measure  a  linear 
relationship  between  the  natural  logarithm  of  loss  terms  and  the  inverse  of 
absolute  temperature  for  Na-montmorillonite  samples.  Their  results  are  repro¬ 
duced  in  Figure  27  (a)  for  a  frequency  of  9.8  GHz  and  in  Figure  27  (b)  for  a 
frequency  of  0.1  MHz.  The  activation  energy  associated  with  the  low- 
frequency  measurements  was  about  12  kcal/mole,  while  that  associated  with 
the  high-frequency  measurements  was  about  6  kcal/mole  as  long  as  the  temper¬ 
ature  exceeded  -52  °C.  The  high-frequency  results  are  consistent  with  the  free 
water  dielectric  relaxation  mechanism  discussed  earlier  and  attributed  to  the 
breaking  of  hydrogen  bonds.  However,  the  low-frequency  measurements  indi¬ 
cate  another  mechanism  for  energy  dissipation,  most  likely  free  charge  carriers 
coupled  with  Maxwell-Wagner  effects. 

Campbell  (1988)  used  the  same  approach  as  did  Hoekstra  and  Doyle  to 
estimate  the  activation  energy  in  soils  above  freezing  temperatures.  Data  he 
measured  at  1  MHz  are  reproduced  in  Figure  28.  Campbell  calculated  an 
activation  energy  of  about  7  kcal/mole  and  argued  that  the  dominant  loss 
mechanism  at  this  frequency  was  the  Maxwell-Wagner  effect. 
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4  Models  for  Soil  Electrical 
Behavior 


Moist  soils  are  a  heterogeneous  mixture  of  mineral  particles,  water,  air,  and 
possibly  any  number  of  organic  substances.  The  mineral  particles  cover  a 
broad  range  of  sizes,  varying  from  as  small  as  a  few  micron.s  to  several  cen¬ 
timeters.  Water  is  not  pure  and,  as  discussed  in  the  previous  section,  in  com¬ 
bination  with  the  structural  features  of  the  mineral  particles,  can  show 
anomalous  dielectric  behavior  at  radio  fret)uencies  unlike  that  of  pure  water. 
Nevertheless,  the  only  practical  approach  to  predicting  soil  response  to  either 
physical  or  electrical  driving  forces  is  that  a  real  soil  must  be  modeled  as  a 
continuum  with  homogeneous  properties. 

As  with  the  review  of  related  electrical  property  measurements  in  the  previ¬ 
ous  chapter,  a  review  of  electrical  property  models  for  soils  and  water  would 
be  remiss  without  a  summary  table  (Table  2)  of  significant  contributions. 

Many  of  the  entries  in  the  table  will  be  expounded  on  in  succeeding  sections. 


Mixing  Models 

The  most  common  approach  to  modeling  the  dielectric  response  of  moist 
soils  that  attempts  to  draw  upon  a  combination  of  the  physical  structure  of  the 
soil  and  expc;1mental  observations  of  the  electrical  response  of  laboratory 
samples  is  that  of  mixnjrc  theory  in  which  the  effective  dielectric  coastant  for 
the  material  takes  the  functional  form 


W  =  [e,  W . 0] 


(33) 
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where 


tjf)  =  frequency-dependent  dielectric  constant  of  the  mixture 

t,(f)  =  frequency-dependent  dielectric  constant  of  the  ith  component  of  the 
mixnire 

In  general,  the  mixture  dielectric  should  lie  somewhere  between  that  of  a  paral¬ 
lel  plate  capacitor  that  is  filled  with  a  mix  of  fibers  extending  from  one  plate 
to  the  other  and  that  of  a  capacitor  filled  with  sheets  of  material  whose  inter¬ 
faces  are  parallel  to  the  plates  (Hasted  1973).  The  former  is  called  parallel 
mixing,  and  the  latter  is  called  series  mixing.  For  example,  if  one  took  the 
soil  to  consist  of  air,  water,  and  some  homogeneous  soil  mineral,  the  two 
mixing  models  would  take  the  form 

PARALLEL  MIXING 

E,  =  H'A  -  .  W,t,  (34) 


SERIES  MIXING 


(35) 


where  the  subscripts  are  self-explanatory  and  the  W,  represent  weighting  fac¬ 
tors.  The  weighting  factors  are  most  commonly  taken  to  be  the  volume  frac¬ 
tion  occupied  by  that  substance. 

One  reference  to  an  early  generalization  of  the  mixture  model  for  three 
components  was  given  by  Ansoult,  deBacker,  and  DeQercq  (1984)  as 

In  this  case,  /t  =  +1  represents  parallel  mixing,  and  k  =  -\  represents  series 
mixing. 

From  an  even  more  physical  basis,  many  researchers  have  developed  dielec¬ 
tric  constant  mixture  models  that  consider  the  effects  of  local  electric  fields  on 
the  induced  internal  fields  of  specially  shaped  homogeneous  dielectric  parti¬ 
cles.  Reynolds  and  Hough  (1957)  argued  that  a  general  formulation  for  a 
three  component  system  could  be  written  as  either 
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(37) 


or 


(£.-£.)5/.  *  (£«-<^)V2  =  0 


(38) 


where  the  5|  are  the  volume  fractions  of  material  "i"  and  the  /  are  the  local 
field  ratios;  that  is,  the  ratio  of  the  held  within  the  dielectric  component  to  the 
average  field  in  the  sample.  From  Stratton  (1941),  a  homogeneous  electric 
field  induces  a  field  inside  a  spheroid  in  such  a  way  that 


/.  =  E 


i< 


co^a- 


rn 

e 

1 

-1 

1*^; 

(39) 


where 

a,  =  angles  between  the  field  and  the  axes  of  the  spheroid 

Aj »  depolarization  factors  related  to  the  shape  of  the  spheroid 

e*  s  dielectric  constant  of  the  homogeneous  material  into  which  the 
spheroid  is  inserted 

deLoor  (1968)  published  the  same  result,  but  less  general  in  the  sense  that  he 
assumed  a  random  orientation  of  sphenoids  for  which  cos^a^  =  1/3. 

One  of  the  more  well-known  mixture  formulae  is  that  associated  with 
Boettcher  which,  for  a  two-component  mixture,  can  be  written  as 

=  S  (40 

K  ‘(£^26,) 


which  comes  from  Equations  37  and  39  by  setting  =  1/3  (spheres)  and  £*  = 
e„.  In  this  formula,  t,  is  the  dielectric  constant  of  the  particles  and  is  the 
dielectric  constant  for  the  host  medium. 

Another  popular  formula  comes  from  Equations  38  and  39  and  e*  =  £3, 

Aj  =  1/3  (spheres).  This  yields  Rayleigh’s  formula 
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^  s 


(41) 


Bruggeman  took  the  differential  form  of  Rayleigh's  formula  and  integrated 
it  between  the  limits  of  and  e„  to  arrive  at  the  formula 


=  e. 


(42) 


Looyenga  (1965)  converted  Boettcher’s  formula  into  a  differential  equation 
on  volume  fraction  8j  and  integrated  between  the  limits  of  0  and  1  to  give  the 
symmetrical  form  (similar  to  Equation  36) 


e  j  =  6.£,^  - 


(43) 


Birchak  ei  al.  (1974)  utilized  a  symmetrical  form  much  like  Equation  43 
except  that  the  powers  were  taken  to  be  1/2  instead  of  1/3. 


(44) 


One  positive  feature  of  the  mixture  models  is  that  they  all  come  from  the 
same  origin;  namely,  consideration  of  the  average  electric  field  and  electric 
displacement  in  a  mixture  of  otherwise  homogeneous  materials  in  which  the 
variables  are  physically  measurable  quantities  of  volume  fractions  and  the 
dielectric  constants  of  the  components.  Thus  there  is  much  physics  involved 
in  developing  these  formulae.  The  approximations  involved  in  obtaining  dif¬ 
ferent  formulae  include  the  particle  shapes  and  orientation  distribution  and 
what  value  gets  assigned  to  the  host  medium  dielectric  constant  when  comput¬ 
ing  tire  field  ratios. 

All  of  these  mixture  formulae,  as  well  as  many  others  (Reynolds  and 
Hough  1957)  have  been  applied  to  particular  data  sets  with  satisfactory  results. 
However,  saying  that  one  formula  applies  to  all  situations  and  materials  would 
be  impossible. 

The  converse  is  also  true;  i.e.,  one  cannot  expect  all  of  tlicse  formulae  to  fit 
a  limited  set  o;  data.  Wang  and  Smugge  (1980)  exercised  several 
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iwo-component  mixture  models  (probably  resulting  in  slight  overestimates  at 
low-moisture  contents  due  to  ignoring  air  voids)  to  calculate  the  real  part  of 
the  complex  dielectric  constant  as  a  function  of  moisture  content  and  to  com¬ 
pare  the  span  of  model  predictions  with  real  data  collected  at  1.4  GHz.  Fig¬ 
ure  29  shows  the  results  of  those  calculations.  Curve  3  is  the  bounding  simple 
parallel  mixing  formula  (Equation  34);  curve  1  is  Rayleigh’s  formula  (Equa¬ 
tion  41);  curve  2  is  Boettcher’s  formula  (Equation  40);  and  curve  4.  which 
gives  the  best  fit,  is  Birchak’s  model  (Equation  44).  Looyenga’s  results  (Equa¬ 
tion  43)  have  also  been  added  as  “L"  symbols. 

Campbell  (1988)  has  demonstrated  that  the  dielectric  response  of  real 
materials  is  bounded  by  the  series  and  parallel  models.  It  is  not  suiprising  that 
1-ooyenga’s  model  for  spherical  particles  is  a  good  fit  to  much  of  the  data 
(Hasted  1973),  as  it  clearly  lies  between  these  limits.  Dobson  et  al.  (1985)  has 
carried  mixture  models  to  the  four-component  level  by  introducing  different 
parameters  for  soil  particles,  air,  free  water,  and  water  lightly  bound  to  the  soil 
particles. 


Equivalent  Circuits 

Thinking  of  mixtures  in  terms  of  their  electrical  analogs  is  often  advanta¬ 
geous,  particularly  v/hen  one  wants  to  consider  different  electrical  paths 
through  the  medium  and  some  semiempirical  way  to  account  for  the  weighting 
of  those  paths.  The  following  paragraphs  take  a  careful  look  at  equivalent 
circuits  for  both  individual  components  of  mixtures  and  the  mixmres  them¬ 
selves.  They  range  from  quite  simple  models  to  quite  complex  circuits.  How¬ 
ever,  with  today’s  readily  available  compuuiional  power,  exercising  these 
models  is  very  straightforward. 


Homogeneous  materials 

The  complex  dielectric  constant  and  conductivity  of  a  linear  material  relate 
an  applied  electric  field  to  current  flow  within  the  material  through  Ampere’s 
law  (Equation  10).  Current  flow  in  the  material  can  be  modeled  by  a  two- 
terminal  network  of  lumped  circuit  elements.  Furthermore,  because  real  mate¬ 
rials  do  exhibit  both  energy  storage  (capacitive)  and  energy-loss  (free  charge 
movement  and  dielectric  relaxation)  mechanisms,  imagining  that  simple 
equivalent  circuits  will  have  to  include  both  capacitive  and  resistive  elements 
is  natural.  The  following  paragraphs  describe  models  of  the  shunt  response  of 
two-terminal  equivalent  circuits  that  represent  the  flow  of  current  (both  con¬ 
duction  and  displacement)  between  the  inner  and  outer  conducting  surfaces  of 
a  coaxial  device  and  through  the  dielectric  mixture  being  measured.  The  series 
response,  or  current  flow  along  the  inner  and  outer  conducting  surfaces,  which 
would  be  represented  by  series  inductive  and  resistive  elements,  is  not  mod¬ 
eled. 
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VOLUMETMIC  WATER  COMTCNT 


Figure  29.  Several  mixing  formulas  ploned  against  real  data  (from  Wang  and 
Smupge  (1980)) 

As  examples  of  simple  equivalent  circuits,  one  should  consider  the  two 
circuits  shown  in  Figure  30,  one  having  constant  value  elements  in  senes  and 
the  other  with  elements  in  parallel.  Drawing  upon  the  very  useful  concept  of 
equivalent  .apacitances  for  each  of  the  elements  in  these  and  the  following 
circuitsCQ  =  i/o)/?;  Q  =  -//(O^L).  the  frequency-dependent  dielectric  response 
of  the  circuits  in  Figure  30  can  be  shown  to  be 
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Figure  31 .  Debye  equivalent  circuit  and  its  response 
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1  ♦{mftC)  J 


when  C„  represents  the  vacuum  capacitance. 
SIMPLE  PARALLEL  MODEL 


E'  =  £  (46) 


E  =  _ 

CflO)/? 


The  physical  interpretation  of  the  parameters  contained  in  these  equations  is 
that  the  capacitor  accounts  for  the  dielectric  polarisation  of  the  material  and 
the  resistor  is  representative  of  ionic  conduction  or  other  losses  that  decrease 
with  increasing  frequency.  Both  parameters  arc  adjustable  to  provide  the  best 
fit  of  the  amplitude  and  phase  response  of  the  material  as  it  is  measured  over  a 
broad  range  of  frequencies. 

However,  neither  of  these  circuits  adequately  model  the  high-frequency 
behavior  of  water  with  regard  to  the  real  part  of  the  dielectric  constant.  The 
parallel  model,  when  placed  in  series  with  another  parallel  model,  has  been 
used  to  model  the  Maxwell- Wagner  losses  in  saturated  soils  (Smith  1971; 
Bidadi,  Schroeder,  and  Pinnavaia  1988). 

As  reported  in  an  earlier  section.  Cole  and  Cole  (1941)  used  equivalent 
circuits  to  describe  both  the  ideal  Debye  relaxation  model  and  their  own  inter¬ 
pretation  of  real  data.  To  help  make  this  section  of  the  study  self-contained, 
their  circuits  arc  reproduced  in  Figures  31  and  32. 

The  corresponding  governing  equations  for  the  Debye  and  Cole-Cole  mod¬ 
els  arc 

DEBYE  MODEL 


e' 


qWa)/?C,)’' 


(47) 
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Figure  32.  Cole-Cole  equivalent  circui*  and  its  response 
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COLE-COLE  MODEL 
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where  l/T,  is  the  radial  frequency  at  the  peak  loss. 


Approximating  the  dielectric  response  of  impure  water  through  the  use  of 
equivalent  circuits  is  even  possible,  von  Hippel  (1954)  used  the  circuit  shown 
on  Figure  33  as  an  approximation  to  the  frequency  response  of  water  exhibit¬ 
ing  low-frequency  ionic  conductivity.  The  frequency  response  of  the  dielectric 
constant  in  von  Hippel’s  water  model  is  governed  ty  the  following  set  of 
equations. 


VON  HIPPEL’S  WATER  MODEL  (with  ionic  conductivity) 


1 - 1 

*  wqj 


1 

Co(u/?;[i  *  [i^ff 
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Figure  33.  Equivalent  circuit  for  impure  water  (after  von  Hippel  (1954)) 

In  fact,  there  is  no  reason  not  to  include  another  parallel  capacitor  to  von 
Hippel's  water  model  to  account  for  the  high-frequency  permittivity,  ending  up 
with  a  four-element  model  shown  in  Figure  34  for  which  the  governing  equa- 
tioas  arc 


t  -  _ -  +  - 


c, 


1  ^ 


(50) 


1^1  -t 

r 


c,(uR.[i 


Mixtures 

it  certainly  seems  reasonable  to  use  combuiations  of  homogeneous  m-atenal 
equivalcni  circuits  to  model  the  frequency  dependent  dielectric  response  of 
mixtures  surh  as  moist  soils  In  fact,  such  models  would  be  particularly  useful 
for  studying  the  effects  of  water  content  in  nonsaiuratcd  soils  if  there  could  be 
a  sensible  weighting  of  the  equivalent  circuit  elements  to  account  for  varying 
moisture  contents. 
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Figure  34.  Water  model  including  DC  conductivity  and  optical  permittivity 


Campbell  (1988)  demonstrated  quite  clearly  that  if  one  considers  the  behav¬ 
ior  of  the  real  part  of  the  complex  dielectric  constant  as  a  function  of  water 
content  for  soils,  it  falls  between  the  limits  for  a  thrce-componcni  system  (soil, 
water,  and  air)  created  by  assuming,  from  one  viewpoint,  that  the  three  compo¬ 
nents  can  be  reprcsented  by  equivalent  capacitances  in  series  and,  from  another 
viewpoint,  by  equivalent  capacitances  in  parallel. 

Bidadi,  Schroedcr,  and  Pinnavaia  (1988)  extended  this  concept  a  bit  further 
to  successfully  model  the  complex  dielectric  response  of  clay  films  as  a 
Maxwcll-Wagner  mechanism  using  the  equivalent  circuit  shown  in  Figure  35. 
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The  upper  frequency  limit  of  their  study  was  only  about  KK)  MHz.  Their 
reasoning  for  choosing  such  an  arrangement  of  elements  is  the  following. 
Tlicy  wanted  to  model  the  physical  situation  of  easy  conducting  paths  through 
the  clay  interstices  that  could  be  interrupted  by  thin  insulating  barriers  (clay 
platelets,  for  example),  this  being  their  interpretation  of  the  Maxwell-Wagner 
effect.  Small  values  of  R,  and  large  values  of  C,  allow  for  ea.sy  current  flow, 
while  large  values  of  R2  and  small  values  of  C2  account  for  the  barriers.  A 
similar  model  was  proposed  by  Eicke  et  al.  (1986)  in  which  the  barriers  (actu¬ 
ally,  the  solvent)  were  taken  as  lossless  (/?2  =  0). 
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Figure  35.  Equivalent  circuit  for  clay  film  studies  (adapted  from  Bidadi, 
Schroeder,  and  Pinnavaia  (1988)) 


Backing  up  further  in  time,  Sachs  and  Spiegjer  (1964)  and,  later.  Smith 
(1971)  applied  even  more  complicated  equivalent  circuits  to  the  study  of  the 
complex  dielectric  response  of  saturated  materials.  Theirs  was  a  very  physi¬ 
cally  satisfying  approach  to  the  problem  in  which  they  assumed  a  combination 
of  equivalent  circuits  in  series  representing  a  discontinuous  electrical  path 
through  the  soil-water  system  (Bidadi ’s  reasoning  for  interrupted  current  flow) 
and  of  equivalent  circuit  elements  in  parallel  representing  continuous  paths  for 
effective  current  flow.  The  resulting  combined  circuit  is  shown  in  Figure  36 
and  was  applied  to  data  that  did  not  exceed  100  MHz. 

The  geometrical  parameters  a,  b,  c,  and  d  were  empirically  established  by 
Sachs  and  Spiegler  for  their  work  by  obtaining  a  best  fit  to  low-frequency 
conductivity  data.  One  has  to  believe,  however,  that  these  parameters  must  be 
related  to  physical  quantities.  For  example,  it  seems  likely  that  d  could  be 
associated  with  porosity.  For  saturated  media,  the  higher  the  porosity,  the  less 
chance  exists  for  series-like  interruptions  in  effective  current  flow.  Likewise,  b 
should  be  related  to  dry  density.  The  greater  the  sample  dry  density. 
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Figure  36.  Equivalent  circuit  model  for  saturated  media  (adopted  from  Sachs 
and  Spiegler  (1964)  and  Smith  (1971)) 


the  more  opportunity  there  is  for  conductive  paths  through  the  soil  matrix. 

And  finally,  t  should  somehow  relate  to  free  water  content  in  the  soil  sample. 
The  greater  the  amount  of  ftee  water,  the  greater  would  be  the  opportunity  for 
current  to  flow  through,  unimpeded  by  the  solid  matrix.  The  governing  equa¬ 
tions  for  the  circuit  representation  in  Figure  36  are 

THE  THREE-PATH  EQUIVALENT  CIRCUIT 

,, . (5]) 
C,(o,’  •  D,>)  C. 


„  __  {N,D,  -  ^  ^ 


where 

N.  =  (a^RR  CC  -  1 

/  S  w  i  w 
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yVj  =  (0  (/?,C,  +  R^cjj 


D,  =  (s9RR^  [{l-d)C^  +  dC^ 

Dj  =  0)  [d«,  +  (l  -4/?J 

Including  another  set  of  equivalent  circuits  to  those  of  Figure  36  to  account 
for  air  in  nonsaturated  soils  certainly  seems  reasonable.  Perhaps  even  more 
simply,  allow  a  +  b  +  c  <  1  and  1  -  d  to  be  replaced  by  x  -  d  where  x  <  1 . 

Furthermore,  arranging  a  limited  number  of  electrical  circuit  elements  in 
such  a  way  as  to  approximate  almost  any  dielectric  response  in  a  single  homo¬ 
geneous  material  is  undoubtedly  possible.  What  is  really  needed  is  an  equiv¬ 
alent  circuit  for  which  the  elements  not  only  allow  a  good  reproduction  of 
spectral  response  over  a  broad  range  of  frequencies,  but  also  have  sensible 
physical  meaning  or  whose  values  can  be  measured  in  the  laboratory. 

Percolation  Transition  or  Long-Range 
Connectivity 

As  pointed  out  in  Chapter  3,  experimental  data  on  moisture  content  varia¬ 
tions  in  soils  shows  that  the  real  part  of  the  dielectric  constant  at  low 
frequencies  increases  more  rapidly  with  increasing  moisture  at  higher  levels  of 
moisture  content  than  it  does  at  lower  moisture  levels.  These  data  can  be 
viewed  as  bilinear  in  moisture  content,  where  below  some  critical  value  of 
moisture  content,  the  rate  of  increase  in  dielectric  constant  is  nearly  constant. 
Above  that  critical  value,  the  rate  is  again  nearly  constant  but  somewiiat 
higher. 

Wang  and  Smugge  (1980)  argued  that  this  transition  point  was  a  function 
of  soil  texture;  i.e.,  a  function  of  the  amount  of  sand  and  clay.  They  also 
related  the  wilting  point  of  the  soil  (the  volumetric  water  content  under  a  pres¬ 
sure  of  15  bars  at  standard  temperature)  to  it.;  texture  and,  in  so  doing,  were 
able  to  show  a  strong  correlation  between  the  transition  point  and  wilting 
point  In  fact,  Wang  and  .Smugge  went  on  to  develop  these  observed  soil 
texture  effects  .nto  empirical  relationships  for  the  electrical  properties  of  soils. 

Campbell  (1988)  viewed  this  phenomenon  by  electrical  analogy  as  a  transi¬ 
tion  between  series  mixing  (Equation  35)  at  lower  water  contents  and  parallel 
mixing  (Equation  34)  at  moisture  levels  above  the  critical  value.  He  borrowed 
from  the  percolation  theory  the  concept  of  a  percolation  transition  in  moisture 
at  which  level  one  or  more  continuous,  low  resistance  electrical  paths  through 
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the  porous  media  suddenly  appear.  He  referred  to  this  transition  as  the  appear¬ 
ance  of  “long-range  connectivity." 

Campbell  attempted  to  create  a  computerized  model  of  moist  soil  passing 
through  the  transition  point  using  the  grid  shown  in  Figure  37.  TTie  grid  was 
initially  set  up  with  a  random  distribution  of  soil  and  air  squares  to  provide  a 
given  porosity.  The  air  squares  were  then  randomly  replaced  with  water 
squares  as  water  content  was  allowed  to  increase.  A  dielectric  constant  was 
then  calculated  for  each  column  of  soil,  air,  and  water  squares  extending  from 
one  capacitor  plate  to  the  other  by  the  series  model  (Equation  35),  and  finally, 
an  effective  dielectric  constant  for  that  moisture  content  was  calculated  as  the 
average  value  for  all  of  the  columns. 


Figure  37.  Campbell's  (1988)  percolation  model  grid  for  moist  soil 

Wliile  this  model  simulated  the  inciieasing  rate  of  change  in  teal  dielectric 
constant  with  increasing  water  content,  it  did  not,  in  fact,  simulate  long-range 
connectivity  because  it  was  only  a  one-dimensional  model.  Unless  the  dry  soil 
grid  contained  several  columns  of  only  air  spaces,  simulating  long-range  con¬ 
nectivity  was  not  possible.  Campbell  also  pointed  out  that  another  shortcom¬ 
ing  of  this  model  was  its  inability  to  scale  properly.  Doubling  capacitor 
dimensions  gave  a  different  result. 


Fractal  Models  of  Electrical  Behavior 

Modeling  of  electrodes 

In  recent  years,  fractals  (sec  Appendix  C  for  deFinitions)  have  been  utilized 
to  study  the  electrical  behavior  of  fluid-saturated  porous  electrodes.  The  fun¬ 
damental  thought  behind  the  use  of  fractal  concepts  is  that  the  number  of  ele¬ 
ments  used  to  construct  a  fractal  geometry  scales  as  a  power  law  of  the  cle¬ 
ment  length  or  relative  size  fMandcIbrot  1983)  It  seems  to  have  become 
fashionable  to  also  talk  about  a  fractal  geometry  as  one  for  which  “physical 
quantities  scale  as  a  power  law  of  the  length  scale. ...on  which  they  arc 
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measured”  (Brouers  and  Ramsamugh  1988).  For  example,  if  there  is  a  region 
of  the  electromagnetic  spectrum  where  conductivity  of  a  sample  of  material  is 
observed  to  vary  as  a  power  law  in  frequency  (hence,  wavelength),  then  that 
sample  is  thought  to  possess  a  fractal  geometry.  In  fact,  Brouers  and 
Ramsamugh  (1988)  combined  low-frequency  permittivity  and  conductivity  data 
for  brine-saturated  alumina  ceramics  with  some  statistical  models  to  calculate 
fractal  dimensions  of  2.65  ±  0.10  for  five  different  ceramic  samples,  indepen¬ 
dent  of  sample  porosity. 

As  part  of  a  study  on  the  frequency  response  of  porous  battery  electrodes, 
Sapoval,  Chazalviel,  and  Peyriere  (1988)  also  observed  that  the  impedance  of 
rough  or  porous  electrodes  often  goes  like  a  power  law  in  frequency  (constant 
phase  angle  behavior)  and  is,  in  fact,  related  to  the  fractal  dimension  of  the 
electrode  surface  in  some  cases.  The  model  they  used  to  associate  fractal 
dimension  to  electrode  structure  was  that  called  a  “finite  modified  Sierpinski 
electrode"  as  shown  in  Figure  38.  The  Sierpinski  electrode  is  something  like  a 
one-dimensional  Menger  sponge  that  is  discussed  in  Appendix  C.  Sapoval. 
Chazalviel,  and  Peyriere  obtained  an  expression  for  the  admittance  of  a  single 
pore  within  the  electrode  and  used  the  fact  that  the  total  pore  admittance  is  the 
sum  of  the  admittances  of  all  the  pores  to  analytically  study  the  effect  of  frac¬ 
tal  dimensions  on  the  electrode  response. 

Dissado  and  Hill  (1989)  used  the  same  model  to  simulate  analytical 
response  functions  and  to  attempt  to  prove  that  there  is  one  form  of  self- 
similarity  at  high  frequencies  (within  individual  fractal  pores)  and  another  at 
low  frequencies  (sum  of  all  pores). 


Fractal  pore-fllling  model 

Campbell  (1988)  observed  that,  up  to  the  time  of  his  publication,  fractal 
concepts  had  been  applied  to  the  bound  water  and  saturated  regimes  for  porous 
media,  but  not  to  the  broad  regime  of  moisture  contents  below  saturation.  He 
was  searching  for  a  way  to  model  the  dielectric  response  of  moist  soils  at 
moisture  contents  less  than  saturation  and  decided  that  it  might  be  possible  to 
do  so  using  fractal  techniques.  His  assumptions  were  as  follows: 

a.  Water  within  pore  spaces  accumulates  in  uniformly  sized  fractal 
clusters. 

The  pore  spaces  are  all  of  equal  characteristic  size  /,  which  is  equal  to 
the  average  particle  size  of  the  soil. 

c.  The  dielectric  constant  of  the  soil-water-air  mixture  increases  with 
increasing  moisture  content  in  a  series-mixing  fashion  until  the  size  of 
the  water  clusters  reaches  the  size  limitation  on  pores  at  which  time 
long-range  water  connectivity  in  the  soil  mixture  is  established;  i.e.,  the 
percolation  threshold  is  reached. 
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Figure  38.  The  “finite  modified  Sierpinski  electrode"  (adopted  from  Sapoval  et 
al.  (1988)) 

d.  Once  long-range  connectivity  is  established,  additional  moisture  enters 
in  a  nonfractal  manner  and  the  dielectric  constant  rises  in  a  parallel¬ 
mixing  fashion. 

e.  The  soil  is  neither  a  clean  gravel  (which  precludes  water  duster  span¬ 
ning  of  the  pone  spaces)  rwr  a  pure  clay  (dominated  by  bound  water 
effects). 

The  key  assumptions  made  by  Campbell  that  allowed  a  rather  simple 
expression  between  critical  water  content  (the  moisture  content  at  which  long- 
range  connectivity  is  established)  and  soil  porosity  to  be  developed  were 
assumptions  a  and  b.  From  a,  Campbell  wrote  that  the  volume  of  water  con¬ 
tained  in  a  cluster  of  size  r  and  fractal  dimension  D  was 

[cluster  water  volume)  o 

(52) 


Thus  the  volumetric  water  content  also  goes  as  r'’.  Mandelbrot,  himself,  found 
no  argument  against  the  mass  density  relationship  that  stated  that  mass  within 
a  fractal  cluster  of  size  R  went  as  (Mandelbrot  1983)  and  if  one  accepts  that 
individual  water  molecules  all  have  the  same  size  and  mass,  then  Equation  52 
is  believable. 
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(53) 


{duster  mass)  o. 


At  such  time  as  the  clusters  just  span  the  pores  of  size  /  to  produce  long- 
range  connectivity,  the  volumetric  water  content  can  be  written  as 


e 


crihco/ 


a  I  ° 


(54) 


Campbell  took  the  soil  sample  porosity  0  to  be  defined  in  three-space  as 


({>  a 


(55) 


He  finally  concluded  from  the  ratio  of  these  last  two  quantities  that 


A  a  (56) 

<t) 


and  that  a  log-log  plot  of  versus  /  should  be  a  straight  line  with  slope 
equal  to  D-3.  If  there  is  no  difference  in  how  the  fractal  water  clusters  grow 
in  different  soils,  then  D  is  the  same  for  any  type  of  soil  that  fits  within  the 
bounds  presented  by  assumption  e.  In  fact,  Campbell  plotted  critical  water 
content  data  for  five  soils  at  a  frequency  of  50  MHz  (shown  in  Figure  39)  and 
found  a  fractal  dimension  of  2.4  that  is  similar  to  others'  calculations  of  about 
2.5  for  aggregate  clustering.  He  made  no  attempt,  nor  is  there  sufficient  data, 
to  draw  any  conclusions  about  the  effects  of  soil  texture  or  signal  frequency  on 
critical  water  content. 

Campbell’s  fractal  pore-filling  model  is  very  simplistic,  but  in  his  own 
words,  “..the  similar  dielecuic  behavior  of  soils  with  huge  ranges  in  particle 
size  and  shape  of  particles  argues  against  a  complicated  model  which  considers 
many  different  soil  parameters..’’  (Campbell  1988).  He  went  on  to  admit  that 
the  weakest  point  in  his  model  was,  in  fact,  assumption  b\  i.e.,  that  the  charac¬ 
teristic  pore  size  equals  the  median  panicle  size  of  the  soil. 
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SCALED  CRITICAL  WATER  COMTEWT  VS.  PORE  SIZE 


Figure  39.  Scaled  critical  water  content  versus  pore  size  (trom  Campbell 
(1988)) 
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Appendix  A 

Skin  Depth  Calculations 


Obviously,  skin  depth  effects  preclude  the  utility  of  practical  electromagne¬ 
tic  devices  to  penetrate  to  distances  comparable  with  acoustic  devices,  but 
nevertheless,  low-frequency  devices  have  some  utility  (Berlin  et  al.  1986;  Farr 
et  al.  1986;  Schaber  et  al.  1986).‘  As  a  means  of  quantifying  the  penetration 
of  electromagneuc  waves  into  earth  media  in  terms  of  the  complex  die’ectric 
constant  of  the  media,  one  should  consider  the  following  argumenus.  The  one¬ 
dimensional  propagation  of  an  electromagnetic  wave  into  some  medium  can  be 
described  by  an  amplitude  function  such  as 


g  >{**■«*) 


(Al) 


where 

{  =  symbol  designating  the  imaginary  quantity 
k  =  complex  wave  number 
X  =  distance 
0)  =  radial  frequency 
t  =  time 

Furthermore,  let 

k  =  P+i  ^  (A2) 


'  Refrrcnces  ciied  in  this  appendix  arc  located  at  the  end  of  the  main  text. 
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where 


A2 


P  =  phase  constant  and  governs  how  a  lossless  medium  propagates  a  wave 

a  =  the  attenuation  factor  and  controls  how  the  amplitude  of  tiie  wave 
decreases  with  distance  traveled  through  the  medium. 

Then  the  amplitude  of  the  wave  is  as  follows; 


i.e.,  a  traveling  wave  combined  with  an  attenuation  factor. 

If  skin  depth  is  defined  as  that  distance  over  which  tlie  wave  amplitude 
decreases  by  a  factor  of  He,  then  the  skin  depth  5  is 


6 


(A4) 


But  the  wave  number  can  also  be  written 

^  ^  0)  _  oW 
V  c 


(A5) 


where 

V  =  phase  velocity  in  the  medium 

N  =  complex  index  of  refraction 
=  n'  +  in" 

c  =  the  speed  of  light 
Tlius 

6=1  =  —  (A6) 

a  ton" 
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or  the  skin  depth  is  inversely  proportional  to  the  imaginary  part  of  the  complex 
index  of  refraction. 

However,  this  study  focuses  on  the  complex  dielectric  constant  e  as  the 
measure  of  soil  electrical  characteristics. 


e  =  e'+ie"  s 


(A7) 


Following  several  substitutions  and  simplifications,  one  can  show  that 


6 


e'WHe"/e')  -1 


(A8) 


where  is  the  free-space  wavelength.  A  visualization  of  this  equabon  for  the 
range  of  values  of  permittivity  e'  and  loss  tangent  e'Ve'  normally  found  in  soils 
is  provided  in  Figure  Al.  This  figure  does  not  represent  any  new  science,  as 
the  same  information  has  been  presented  in  other  forms  (von  Hippel  1954).  It 
is.  however,  a  more  condensed  way  of  visualizing  the  skin  depth  that  has  not 
been  observed  elsewhere  in  the  literature.  Figure  Al  is  valid  for  any  interpre- 
taUon  of  loss  tangent:  i.e.,  whether  one  is  speaking  of  a  loss  mechanism  due  to 
actual  charge  migration  or  one  due  to  dielectric  relaxation. 

Another  useful  sldr  depth  relauonship  found  in  the  literature  is  the  nomo¬ 
graph  generated  by  Albrecht  (1966)  and  reproduced  on  Figure  A2.  Albrecht 
first  developed  an  empirical  relabonship  that  associated  ground  conduebvity  to 
gravimetric  moisture  content  and  ground  temperature  at  low  frequencies.  He 
then  related  skin  depth  to  wavelength  and  conductivity  by  the  usual  equation 
for  highly  conducbve  materials  (assumes  that  the  conductivity  phenomenon 
dominates  the  loss  mecharusm  in  soils  at  low  frequencies): 


a  = 


\  TtpCO 


(A9) 


where  the  units  used  are  MKS,  and  the  skin  depth  is  in  meters  and  the  ground 
conduebvity  is  in  units  of  mho/meter.  The  magnetic  permeability  was  taken  to 
be  4icX10'’  Henrys/meter.  The  w'ay  to  use  this  nomograph  is  to  select  a  mois¬ 
ture  content  and  move  on  a  horizontal  line  to  the  appropriate  temperature 
curve.  Then  move  on  a  vertical  line  to  the  appropriate  frequency  curve.  From 
that  point,  a  horizontal  line  to  the  right-hand  scale  reveals  the  desired 
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value  of  skin  depth.  There  is  no  way  to  relate  Albrecht’s  rcsiulis  to  those  of 
Figure  Al,  because  it  begins  with  an  unknown  empirical  relationship  among 
Vmoisture,  temperature,  and  conductivity. 


Figure  A1.  Skin  depth  as  a  function  of  wavelength,  permittivity,  and  loss 
tangent 


Hoekstra  and  Delaney  (1974)  also  published  some  useful  data  on  plane 
wave  attenuation  in  moist  soils  at  higher  frequencies.  These  data  are  shown  in 
Figure  A3.  They  concluded  that  anenuation  was  relatively  independent  of  soil 
fabric  and  that  either  passive  or  active  microwave  sensors  could  obtain  a 
measureable  response  from  only  about  the  first  5  cm  of  ground  below  the 
surface. 
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GRAVIMETRIC  MOISTURE 


Figure  A2.  Skin  depth  nomograph  (from  Albrecht  (1966)) 
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Appendix  B 

Reflection  and  Refraction  at 
Plane  Interfaces 


Remote  sensing  of  the  environment  usually  involves  the  interaction  of  elec¬ 
tromagnetic  energy  waves  with  solid  media  at  some  air-media  interface.  For 
example,  an  active  radar  transmits  waves  toward  a  surface  (usually  at  some 
oblique  angle)  and  measures  what  is  reflected  back  toward  itself  (backscatter) 
or  toward  some  other  receiver  (bistatic  reflection).  Passive  devices  receive 
energy  emitted  by  the  surface  (and  the  atmosphere  within  the  antenna’s  beam- 
width)  due  to  internal  .sources  as  well  as  what  is  reflected  from  natural  sources 
such  as  the  sun,  the  atmosphere,  and  surrounding  terrain. 

An  idealization  of  the  wave  interaction  phenomenon  is  shown  in  Figure  B1 
where  an  incident  plane  wave,  defined  by  its  wave  vector^r ,  gives  rise  to  a 
reflected  wave  k  ^  and  a  refracted  (or  transmitted)  wave  k  j.  (Nearly  every 
textbook  on  electromagnetic  theory  has  a  development  on  reflected  and 
refracted  plane  waves.  What  follows  was  taken  primarily  from  Jackson 
(1975)'  with  a  slight  change  in  notation.  Units  are  cenbmeter-gram-second.) 
The  wave  vectors  are  related  to  the  dielectric  properties  of  both  media  by 

If, I  .  If.!  .  t  .  ^ 


-  (B2) 

|Xj.,  -  -  — —  - - - 


'  References  ciied  in  ihis  appendix  are  located  at  the  end  of  the  mam  lexi. 
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Figure  Bl.  Wave  vectors  at  a  plane  interface 
where 


0)  =  rariiaJ  frequency 

c  =  speed  of  light  in  a  vacuum  (3  x  10*  m/sec) 
p  =  magnetic  pemeability  (hereafter  taken  to  be  unity) 

£  =  complex  dielectric  constant 
N  =  complex  index  of  refraction 

The  plarie  of  incidence  is  taken  to  be  that  plane  formed  by  the  k,  and  n  vec¬ 
tors;  i.e.,  the  plane  of  the  paper. 

Under  these  ideal  conditions  of  plane  waves  at  a  smooth  interface,  one  can 
write  the  following  expressions  for  the  three  sets  of  electromagnetic  energy 
fields. 

INCIDENT 


B2 
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k  iXE 


'1^1 


k,xE, 


k  (0 


REFRACTED  (TRANSMITTED) 


(B4) 


Ej - £76 


1  {kr^-iat) 


(B5) 


■/t  T  X  £ 


kf  c 


(B6) 


REFLECTED 


{kp-x-Uil) 


(B7) 


^R  -  /Mi£i  ^  ^  ~  ^R  ^  ^ R 


(B8) 


By  insisting  that  the  phase  factors  {k  ■  x  terms)  are  all  equal  at  z  =  0,  one  is 
led  to  Snell's  law  and  the  fact  that  6^  =  8;^. 


sin9;  kj 

sinB^-  k 

What  one  desires  to  have  arc  expressions  for  the  amplitudes  of  the  reflected 
and  refracted  waves,  for  it  is  the  wave  amplitude  (or  voltage  amplitude)  that  is 
measured  by  a  receiver  unit.  Amplitude  relationships  can  be  found  by  apply¬ 
ing  the  electromagnetic  boundary  conditions;  namely,  that  the  normal  compo¬ 
nents  of  electric  displacement  Dand  magnetic  induction  B  arc  continuous 
across  the  interface  as  arc  the  tangential  components  of  the  electric  field  L  and 
the  magnetic  field  R.  These  conditions  may  be  written  as  follows; 


_  ^2  (B9) 
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I,  *’  ^  k]  ‘  'I  "  0 


(BIO) 


B4 


[^/  ^  ^ I,  *  ^  ^ K,  ~  ^T  ^  ^T,  *  ”  “  0 


^  -  ^r. 


X  n  =  0 


(BID 

(B12) 


1  (*-,  X  £-,  .  X  £.J  -  1  (.%  X  f ,) 


X  li  =  0 


(B13) 


The  Fresnel  Coefficients 


At  this  point,  two  further  idealizations  are  made  before  finally  writing  down 
the  expressions  for  reflected  and  refracted  wave  amplitudes.  One  is  to  assume 
that  the  incident  electric  field  vector  is  perpendicular  to  the  plane  of  incidence 
(into  the  paper  with  respea  to  Figure  Bl).  and  the  other  is  to  assume  it  is 
parallel  to  the  plane  of  incidence  (lying  within  the  paper  with  respect  to  Fig¬ 
ure  Bl)-  In  radar  parlance,  the  terms  are  “horizontal  polarization”  and  “verti¬ 
cal  polarization,"  respectively. 

Due  to  redundancies  generated  by  Snell’s  law,  one  only  needs  the  third  and 
fourth  boundary  conditions  expressed  in  Equations  B 12  and  B 13  to  arrive  at 
the  final  expressions  for  amplitude  ratios,  which  are  known  as  the  Fresnel 
reflection  coefficients. 

HORIZONTAL  POLARIZATION 


2cos0, 


cos9,  +  —  I  -  sin^0, 

M,e, 


(B14) 
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Reflection  From  Lossless  Media 

The  previous  equations  apply  to  lossy  as  well  as  lossless  media.  For  loss¬ 
less  media,  the  dielectric  constant  and  index  of  refraction  are  real  quantities. 
Refraaed  waves  propagate  without  attenuation,  and  both  refracted  and 
reflected  waves  retain  the  polarization  of  the  incident  wave. 

Of  practical  value  to  the  remote  sensing  arena  is  the  variation  with  incidence 
angle  of  the  reflection  coefficient  in  a  lossless  media.  Because  received  power 
in  a  sensor  goes  as  the  square  of  the  amplitude  of  the  signal  and  is  propor¬ 
tional  to  the  radar  cross  sections  of  the  object  being  illuminated,  scientists  and 
engineers  have  a  simple  model  between  electrical  properties  of  the  media  and 
its  response  to  electromagnetic  energy.  What  is  deduced  from  the  sensor  mea¬ 
surements  and  a  knowledge  of  the  surface  slopes  is  an  apparent  dielectric 
constant  or  apparent  permittivity  of  the  medium. 

Figure  B2  is  a  representation  of  how  the  apparent  permittivity  of  a  medium 
affects  the  reflected  amplitude  ratio  as  a  function  of  incidence  angle  for 
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nonmagn-iic  materials.  The  angle  at  which  the  reflected  wave  is  totally  elimi¬ 
nated  in  the  vertical^jwlarization  case  is  referred  to  as  the  Brewster  angle 
(tanSg  =  N2/N1  =  ^e2/Ei).  For  angles  less  than  the  Brewster  angle,  the  re¬ 
flected  wave  suffers  a  phase  reversal.  The  values  of  dielectric  constant  were 
chosen  to  span  the  range  of  values  from  dry  soil  to  pure  water. 

Reflection  From  Lossy  Media 

Appendix  A  gave  some  indication  of  what  happens  to  electromagnetic 
waves  incident  on  a  lossy  media  (e"  >  O).  The  refracted,  or  transmitted,  wave 
is  attenuated.  Furthermore,  for  oblique  (non-zero)  angles  of  incidence,  sur¬ 
faces  of  constant  amplitude  are  no  longer  parallel  to  the  surfaces  of  constant 
phase,  resulting  in  what  is  referred  to  as  an  inhomogeneous  wave  (Stratton 
1941).  In  fact,  the  index  of  refraction  for  lossy  media  becomes  a  function  of 
the  angle  of  incidence  of  the  incoming  wave.  As  lossiness  goes  up,  the 
refracted  wave  vectors  get  closer  and  closer  to  normal  to  the  interface. 

Similarly,  it  is  possible  to  rationalize  the  second  expressions  in  Equa¬ 
tions  B7  and  B8  to  obtain  ratios  of  reflected  amplitudes  to  incident  amplitudes. 
One  representation  of  these  relationships  may  be  found  in  Stratton  (1941). 

The  effects  of  adding  a  loss  term  to  the  dielectric  constant  are  to  enhance 
reflectance  and  to  reduce  the  destructive  interference  mechanism  that  causes 
the  total  loss  of  reflected  energy  at  the  Brewster  angle.  A  visual  comparison 
of  the  effects  of  loss  factors  in  materials  is  shown  in  Figure  B3. 

It  is  not  inconceivable  that  if  one  had  good  reflectance  data  at  a  single  fre¬ 
quency  over  a  large  range  of  incidence  angles  (enough  to  include  the  minimum 
on  the  vertical  polarization  curve)  and  if  the  wavelength  of  the  sensor  could 
assume  smooth  interface  conditions,  then  one  could  iteratively  obtain  a  first 
approximation  to  the  complex  dielectric  constant  for  the  material  near  the 
surface. 
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Figure  B3.  Reflection  amplitudes  and  phase  shifts  for  nonmagnetic  lossy  mate¬ 
rials 
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Appendix  C 

Fractal  Models  of  Soil 

Structure 


Fractals 

In  recent  years,  some  very  interesting  work  has  been  done  on  the 
conductive  behavior  of  porous  media  using  the  concept  of  self-similar  geome¬ 
tries.  or  fractals,  to  model  the  behavior  of  certain  physical  quantities  of  that 
media.  Mandelbrot,  the  gum  of  fractals,  defines  self-similar  shapes  as  those  in 
which  a  certain  part  of  the  shape  can  be  broken  up  into  N  smaller  parts,  each 
looking  like  the  original,  but  reduced  in  size  by  a  fractional  factor  labeled  r 
(Mandelbrot  1983).'  The  fractal  dimension  D  of  this  self-similar  shape  can  be 
written  as 


D  = 


log  (n) 


V'-l) 


or 


N  =  1 


(C2) 


Mandelbrot  insisted  that  the  exponent  D  be  thought  of  as  a  dimension  because 
it  arose  from  a  method  of  measuring  the  perimeter  of  an  object  having  an 
irregular  boundary. 

Fractals  are  often  used  to  generate  complex,  or  textured,  curves  or  surfaces 
from  very  simple  initiator  geometries.  For  example,  consider  Figure  Cl  in 


'  References  cited  in  this  appendix  are  located  at  the  end  of  the  irnun  text. 


Appendix  C  Fractal  Models  ot  Soil  Structure 


Figure  C1  Fractal  snowflake  of  dimension  1.5  (Mandelbrot  1983) 


which  the  initiator  (or  the  initial  segment  of  the  geometry  that  will  be  fracta- 
lized)  is  a  solid  square  (precisely  stated,  each  side  of  the  square  is  an  initiator), 
and  the  generator  (or  the  desired  geometry)  is  the  broken  line  segment  as 
shown.  Moving  clockwise  from  the  initiator,  the  first  two  applications  of  the 
generator  are  shown  at  full  scale,  while  the  next  two  are  enlarged  to  show  the 
edge  detail.  Figure  C2  shows  another  construction  initiated  with  a  solid  square 
in  which  the  generator  results  in  the  formation  of  both  “lakes”  and  "islanos." 


Soil  Structure 

A  third  example  of  a  fractal  construction  that  is  useful  for  studying  porous 
media,  and  may  be  particularly  useful  for  studying  soil  structure,  is  that  shown 
in  Figure  C3  in  which  the  initiator  is  a  cube  with  square  holes  centered  on 
each  face  and  joining  in  the  center.  The  generator  is  the  same  object  but 
reduced  by  a  factor  of  one-third.  Therefore,  in  the  first  stage,  the  initiator  is 
divided  into  20  smaller  cubes,  resulting  in  a  fractal  dimension  of  2.7268. 
Mandelbrot  referred  to  this  construction  as  a  Monger  sponge.  Friesen  and 
Mikula  (1987)  used  tlie  Monger  sponge  to  model  the  fractal  dimension  of  the 
surface  area  of  coal  in  such  a  way  that  the  slope  of  the  isotherm  (fractional 
volume  change  versus  pressure  for  porosimetry  measurements)  plotted  on  a 
log-log  scale  is  linearly  related  to  the  fractal  dimension.  At  least  one  recent 
study  attempted  to  measure  the  fractal  dimensions  of  sandstone  solid-void 
interfaces  (Katz  and  Thompson  1985).  Measured  values  of  interface  fractal 


C2 
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Figure  C2.  Fractal  islands  and  lakes  of  dimension  1.6131  (from  Mandelbrot  (1983)) 

dimensions  on  five  different  samples  using  electron  microscope  techniques 
varied  from  2,57  to  2.87. 

Work  has  been  done  very  rccenUy  on  modeling  the  microstructure  of  soil 
using  fractals  (Krepfl,  Moore,  and  Lee  1989;  Moore  and  Krepfl  1991).  A 
two-dimensional  representation  of  a  three-dimensional  model  using  hexagonal¬ 
shaped  flakes  as  the  generator  is  shown  in  Figure  C4.  Such  models  may  be 
useful  as  a  representation  of  pure  clays.  This  particular  fractal  stnicture  is 
strongly  reminiscent  of  microphotographs  of  the  pure  kaolinite  mineral 
structure. 
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